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however.  The  second  model  (II)  does  account  for  keel  structure  and  assigns  the 
gaussian  components  of  the  first  model  respectively  to  two  distinct  physical  pro¬ 
cesses:  that  of  local  ice  generation  fprcrcess-S-),  of  S$,  and  that  of  the  inter¬ 
mittent,  large-scale  keel  formation,  [process  K  and  elevation  CjJ.  Two  variants 
of  II  are  also  distinguished  here.  J 

The  empirical  basis  for  these^models  is  the  underice  historgram  of  Lyon  [1], 
[2],  which  can  be  interpreted  to  support  the  concept  of  a  distinct  keel  structure 
(Model  II) , -although  the  precise  existence  of  the  latter  is  subject  to  some  dis¬ 
pute  [3].  On  physical  grounds  it  appears  reasonable  to  treat  the  formation  of 
the  comparatively  small-  and  large-scale  ice  structures  as  separate  physical  pro¬ 
cesses,  each  the  resultant  of  many  effectively  independent  forces,  so  that 
gaussian  elevations  follow  (at  least  approximately),  as  indicated  by  the  Lyon 
data. 

•  In  this  Report  the  analysis  is  limited  to  Model  I  and  calculations  of  back- 
scattering  cross  sections  appropriate  to  intermediate  and  high  frequencies 
.^0(>  0.1  kHz),  i.e.,  Kirchoff-based  methods.  The  treatment  here  is  also  re¬ 
stricted,  for  the  moment,  to  zero  gradient  •(?c*0'j  and  single  ice  bottom  surface 
interactions.  In  most  instances,  the  Rayleigh  numbers  for  both  components  of  the 
underice  surface  are  large,  so  that  theories  based  sojely  on  perturbation  tech¬ 
niques  and/or  Bragg  scatter  are  inadequate.  -  H 

The  construction  of  second-order  characteristic  functions  (c.f.'s)  of  the 
underice  profile  c,  is  accordingly  required.  The  second-order  c.f.  is  modeled 
here  by  a  single  gauss  process,  whose  first-order  moments  are  those  determined 
from  the  two-component,  first-order  gauss  process  fitted  to  the  Lyon  data. 

Because  of  the  large  Rayleigh  numbers  involved,  an  (artificial)  splitting  of  the 
surface  wave  number  spectrum  must  be  made,  to  separate  the  "rough"  (i.e., 
Kirchoff)  components  from  the  Bragg  (i.e.,  "smooth")  components  of  the  scattering 
surfaces.  The  analytic  methods  and  results  used  here  are  direct  modifications 
of  the  author's  recent  studies  of  scattering  from  moving  ocean  wave  surfaces  [4]. 
The  present  treatment  is  directly  applicable  to  intermediate  O(>100  Hz)  and  high 
frequencies  o(>  1-0  kHz). 

In  addition  to  the  analytic  results  for  scattering  cross  sections  developed 
here  for  Model  I,  it  is  found  that  the  simple,  single  surface  version  (Model  la) 
is  inadequate  to  account  for  the  (much)  larger  empirical  (back-)scatter  cross 
sections  observed,  particularly  at  small  grazing  angles  (<j>  =  10°)  and  high  fre¬ 
quencies  (fQ  =  1  kHz).  Accordingly,  a  two-component  version  (Model  lb)  is  pro¬ 
posed,  where  now  an  additional,  independent  surface  projection  rides  on  the 
single  surface  of  Model  la.  RMS  elevations  o(a|<  =  8-15  cm)  and  correlation  dis¬ 
tances  o(Jtc  =  4-7  cm)  for  this  second  surface  component  are  required  for  reason¬ 
able  agreement  with  the  data  cited  [5],  [24].  It  is  emphasized  that,  at  the 
present  stage  of  limited  available  data,  it  is  the  choice  of  parameter  values, 
rather  than  any  specific  analytic  model  form,  which  is  critical  to  the  potential 
applicability  of  the  model  in  question. 

The  results  developed  here  also  indicate  the  need  for:  (1),  profile  eleva¬ 
tion  statistics  with  a  much  finer  resolution  than  presently  available,  say 
o(A0  =  6  m  -  l  cm),  with  particular  emphasis  on  the  range  0(20  cm  -  1  cm); 

(2),  (rms)  slopes  and  slope  statistics;  (3),  directional  data  on  profile  eleva¬ 
tions;  and  (4),  scattering  results  for  Model  II,  both  for  those  cases  where  a 
definite  keel  structure  exists,  and  for  comparison  with  Models  Ia,b  as  generated 
here. 


(Form 


pni  /i7i  \ 


l 


assified 


•I*] 


L 


TD  7375 


TABLE  OF  CONTENTS 


List  of  Illustrations . vii 

List  of  Principal  Symbols . viii 


1. 

2. 


3. 


4. 


Introduction  .  1 

Underice  Surface  Elevation  Models:  I  and  II  .  5 

2.1  Model  I:  An  Unstructured  Two-Component  Underice  Profile  .  5 

2.2  Model  I la :  A  Two-Component  (Two-Dimensional) 

Underice  Profile,  with  Keels  .  10 

2.3  Model  lib:  A  Two-Process  Underice  Profile, 

with  Separate  Keels  and  Bottom  Components  .  13 

2.4  A  Model  for  the  Marginal  Ice  Zone  (MIZ)  (Model  lib)  ...  14 

Scattering  Cross  Sections  for  Model  la  .  14 

3.1  Scattering  Cross  Sections:  Model  la  .  17 

3.2  The  Ref lection^  Tilt-,  and  Shadowing  Factors, 

V  TH-inc’  T . 20 


j\.  The  (Plane-Wave)  Reflection  Coefficient,  r 

The  "Tilt-Factors"  N«j>1ne,  S<°>coh  - 

C..  Shadowing  Factor,  s? . 


3.3  Model  la:  A  Numerical  Example--Backscatter  at 


0oT  =  80°> 


1  kHz 


A.-  Calculations  (0  -  =  80°;  f  =  1  kHz;  <d>, > 2  =  0°-50°)  . 

oi  0  fni  L 

JL  The  Backscatter  Cross  Section,  . 

C,.  Remarks  and  Conclusions . 

Scattering  Cross  Sections  for  Model  lb  . 

4.1  Scattering  Cross  Sections;  Backscatter  Cases,  Model  lb  .  . 

4.2  The  Tilt-Factor  N.^.  „  .  . . 

Lh-inc-back 


4.3  Model  of  the  Small-Scale,  Additive  Surface  Component,  cr 

4.4  Model  lb.  A  Numerical  Example:  Backscatter  @ 

8ot  "  80°»V  a  */2'  fo  =  1  kHz . 

Concluding  Remarks  . 


20 

21 

23 

25 

26 
26 
27 
30 
30 
32 
34 

36 

38 


v 


TD  7375 


Appendix  A.l:  Wave  Number  Spectra  and  Spatial  Correlation 

Functions:  Models  la,  lb .  41 

A. 1-1  One-Dimensional  Correlation  Functions 

and  Spectra  (Model  la) .  41 

A. 1-2  (Isotropic)  Two-Dimensional  Cases 

(Model  la) .  42 

A.  1-3  RMS  Slopes  (Model  la) .  43 

A. 1-4  Covariance  Functions  and  Wave-Number 

Spectra  (Model  lb) .  45 

Appendix  A. 2:  Scattering  Cross  Sections:  Definitions, 

Nomenclature,  and  Conditions  .  50 

A. 2-1  Incoherent  Scattering  Cross  Sections,  3^.  50 

(0lnc 

A. 2-2  Coherent  Scattering  Cross  Sections,  /  .  53 

A. 2-3  The  Reference  Area,  AR£p .  54 

A. 2-4  General  Geometry  in  Detail:  The  Received 

Scattered  Wave .  54 


A. 2-5  Parameters  of  Eq.  (3.6):  R»1  ...  58 

Appendix  A. 3:  Second-Order  Probability  Densities  of  the 

Mixed  or  Ga us s -Compos i te  K,S  Process,  Eq.  (2.1).  .  61 

References .  67 


LIST  OF  ILLUSTRATIONS 


TO  7375 


Figure  Page 

1.1  Historgram  of  underice  profile  5,  from  Lyon's  data  ....  3 

2.1  Roughness  wave  number  (intensity)  spectrum  of  underice 

profile  5,  after  Mellen  et  al . ,  from  FFT  of  Lyon's  data  .  8 

2.2a  Model  II:  Two-dimensional  underice  profiles  (c) 

Eq.  (2.13),  with  distinct  keels  (s^) .  11 

2.2b  Model  I:  Cross  section  of  a  (two-dimensional)  underice 
profile  (?),  without  distinct  keels,  governed  by  Eqs. 

(2.  l)-(2.7) .  11 

3.1  Parameters  of  the  Ice-Bottom  Scattering  Models 

(Models  la,  lb) .  24 

4.1  Backscatter  Cross  Sections  at  0qT  =  80°;  fQ  =  1  kHz  .  .  .  28 

4.2  Sketch  of  Kh(Ar),  Eqs.  (4.12),  (4.12a)  .  37 

A. 2.1  (a)  Backscatter  Geometry  (far-field):  Monostatic 

scattering  (R@T),  vide  Fig.  A. 2. 2;  (b)  Bi static 
Scattering  (Rj*T):  "forward"  or  oblique  scattering 
geometry,  vide  Fig.  A. 2. 2 .  51 

A. 2. 2  Geometry  of  the  transmitter  (@  Oj),  underice^surface  (E), 
receiver  (at  Oq) ,  showing  surface  elevation  X  above 
5  s  0:  S  ,  andKthe  common  region  of  illumination  on  S 
due  to  overlapping  beam  patterns  Q.y,  .  55 


LIST  OF  PRINCIPAL  SYMBOLS 


A  =  beam  parameter 

A ,  =  beam  area  in  illuminated 
1  surface 

°ox- >•  "5 

B  =  beam  parameter 

b,  =  cose„x  +  cose~ 
o  oT  oR 

*  spectrum  parameter 

cT,  Cy  =  speed  of  sound  in  ice, 
water 


kc  =  spectrum  wave  number  parameter 
k  ,  k,  k^  =  wave  numbers 

Lq  =  "bump"  length 

=  correlation  distance 
=  "bump"  length 
A,  =  regions  on  a  plane 
=  source  wavelength 
=  spectrum  splitting  wavelength 

m^,  m^  =  mean  values 


5,6=  unit  functions  on  a 
finite  region 

§_=  (2-dimensional)  delta  function 
5=  delta  function 


F^.  F2  =  characteristic  functions 

fr  =  (carrier)  frequency  of  sound 
source 

fQ  =  spectrum  splitting  frequency 

y  =  density  of  "bumps"  (ice  pro¬ 
jections)  per  unit  surface 
area 


l(2)^  j(4)  _  integrals 


■Ir 


unit 

vectors 


J  =  Bessel  function  of  first  kind 
o 

K^,  =  surface  covariances 

K.,  K2  =  modified  Bessel  functions 
1  c  of  second  kind 


,  ^LH^’  NLh^  =  "t''^t''',^actors', 
rL,  n,  ,  n-  =  unit  normals  to  a 
■~K  ~L  surface  (K.L.S) 


<4>P  *  rms  ice  slope  angle 

pdf  =  probability  density  function 

4>oj  =  horizontal  angle 


R  =  (plane-)  wave  reflection 
coefficient 

R,  R^,  R^  =  rayleigh  number 

jTj ,  r,  Ar,  Ar  =  distances 

P,  Ph,  P^,  Pp,  PK  =  normalized 

covariances 

Pc  =  radius  of  curvature 


S  =  shadowing  function 

3^  =  scattering  cross  section 

aL(x=y)  =  ms  sl°Pes 

Op,  as,  a^,  a^,  =  rms  elevations 

5^  =  amplitude  spectral  density 


vertical  angle 


.  ,  W.,  W„  =  wave  number 

1  c  intensity  spectra 

first-order  pdf 


C,  Cu»  tv*  Cc»  ,  Cu  *  profile 

h  K  S  L  H  elevations, 

components 

Sy*  S|_x »  5Ly  =  profile  slopes 

CQ  =  "bump"  elevation 


TD  7375 


INTERMEDIATE  AND  HIGH-FREQUENCY  ACOUSTIC  BACKSCATTERING 
CROSS  SECTIONS  FOR  WATER-ICE  INTERFACES: 

I.  TWO-COMPONENT  PROFILE  MODELS* 

by 

Dav'd  Middleton** 


I.  Introduction 

The  purpose  of  this  study  is  twofold:  (1),  to  develop  several 
plausible  statistical -physical  models  of  typical  underice  "profiles"t 
of  elevations,  £,  which  can  then  be  used  (2),  to  obtain  acoustic  scat¬ 
tering  intensities  and  spectra.  Here  we  shall  describe  two  general 
models  of  underice  elevations,  and  use  the  simpler  one  for  an  initial 
determination  of  backscatter  cross  sections.  The  extension  of  these 
results  to  the  more  complex  second  model  is  reserved  to  a  subsequent 
study,  for  which  the  present  investigation  is  prototypical.  In  addition 
here,  we  shall  limit  our  analysis  to  zero-gradient  oceans  and  single 
ice-bottom  interactions.  The  emphasis  here  is  on  acoustic  backscatter- 
ing,  but  the  methods  described  are  directly  extendable  to  bi-static 
geometries. 

We  base  both  our  models  on  the  (one-dimensional)  underice  profile 
data  of  Lyon  [1],  [2],  [24].  Figure  1.1  shows  these  data  fitted  by 
two  (first-order)  gaussian  probability  densities  for  £.  No  attempt 
has  been  made  in  this  initial  fitting  to  apply  any  specific  "best-fit" 
techniques:  simple  trial  adjustments  of  means  and  variances,  only,  were 
used  to  obtain  a  visually  reasonable  fit,  which  can  be  improved  upon 
but  without  major  changes  in  the  resulting  parameter  values. 

Our  first  model  (I,  and  specifically  la  here)  postulates  two 


*Work  supported  under  contract  with  the  Naval  Underwater  Systems 
Center,  New  London,  Conn.,  N00140-83-M-PL78,  1983,  and  in  part,  under 
Contract  N00140-84-M-MS42 ,  1984. 

♦♦Contractor:  Physics  and  Applied  Mathematics,  127  E.  91  St., 

New  York,  NY  10128. 

fWe  shall  use  the  term  "profile"  here  interchangeably  with  "surface 
elevation,"  which  is  a  function  of  surface  coordinates  (x',y'),  cf. 

Fig.  (2.2a)ff. ,  e.g.,  r,  =  djr). 
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independent  surface  processes  combining  in  some  fashion  to  produce  a 
single  surface,  c,.  These  two  processes  consist,  respectively,  of  a 
comparatively  small-scale  [m^  =  4.5  m;  a s  =  o(l  m)]  elevation,  in 
turn,  which  is  accompanied  by  a  larger-scale  surface  component,  CK, 

[mK  =  12.5  m;  =  3.0  m] .  As  noted,  these  coexist  in  some  unspeci¬ 
fied  fashion  geometrically,  but  in  such  a  way  that  the  di stributions 
(and  probability  densities)  of  and  add*,  in  the  manner  of  Figure 
1.1,  to  give  the  distributions  (etc.)  of  the  total  profile  C,  cf.  Sec. 

2.1  ff.  This  first  model  (I)  does  not  specify  keel  structure  as  such 
explicitly.  On  the  other  hand,  our  second  model  (II)  does,  by  identi¬ 
fying  the  second  (or  larger-scale)  gaussian  elevation  with  keel  struc¬ 
ture,  on  which  may  "ride"  the  smaller-scale  gaussian  elevation  component 
(Model  lib).  Unlike  Model  I,  the  total  elevation,  or  profile  process 
S,  is  an  (appropriately  weighted)  sum  of  the  independent  processes  £<- 
and  cK,  cf.  Sec.  2.2  ff.  Both  these  models  are  described  analytically 
in  Section  2  following.  Because  there  is  some  doubt  as  to  the  unambig¬ 
uous  geometrical  existence  of  keels  [3],  Model  I  (which  is  the  more 
phenomenological  of  the  two)  may  prove  a  useful  candidate  model.  Decision 
between  them  will  depend  both  on  additional  data  and  on  how  well  the 
corresponding  theoretical  results  agree  with  experiment. 

In  any  case,  we  emphasize  that  both  models  are  potentially  capable 
of  handling  all  frequencies  (above  0 (30-50  Hz))  of  incident  acoustic 
radiation, as  noted  in  Sec.  3  ff.:  neither  is  restricted  to  the  small 
Rayleigh  number  cases  of  Bragg  scatter  alone  [5],  and  each  can  employ 
appropriate  shadowing  and  reflection  functions  [4],  Apart  from  the  fit 
to  the  empirical  data  [cf.  Fig.  1.1],  gaussian  distributions  (and  hence 
gaussian  processes  ultimately)  appear  more  reasonable  than  other  choices. 
In  a  fundamental  way  we  would  say  that  the  process  of  surface  formation 
consists  of  the  superposition  of  many  independent  effects,  so  that  the 
usual  Central  Limit  Theorem  arguments  [6]  apply.  Moreover,  the  presence 


♦Note  that  this  is  not  the  same  thing  as  the  distribution  (and  pdf) 
of  the  process  sum  t,  -  a^s  +  c^‘  ^ec'  ^  See  Appendix  A. 3  ff., 
also. 
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Figure  1.1:  Histogram  of  underice  profile  C,  from  Lyon's  data  [1]; 

Ref.  4  of  [2].  (Solid)  curves:  two  component  gaussian 
fit;  dotted  line:  composite  gauss,  5  =  b.^+b-s,,,  cf. 
(2.7)  et  seq..  1 


of  widely  separated  nonzero  mean  profiles  (elevations),  e.g.,  m<.  = 

4.5  m  vs  =  12.5  m.  Fig.  1.1  (measured  from  the  mean  ocean  surface), 
with  appropriately  small  variances  (o^,o^),  cf.  Fig.  1.1,  ensures  the 


practical  vanishing  of  the  pdf's  w^> 


w^K,  etc.,  in  question  at  c  =  0, 


as  required  (from  the  definition  of  the  profile,  cf.  Fig.  2.1).* 


♦This  latter  requirement,  viz.  c^O,  has  suggested  that  Rayleigh 
processes  for  may  provide  appropriate  models,  but  direct  attempts 

to  construct  Model  I  using  such  processes,  as  well  as  linear  combinations 
of  £e  and  z,v  (Model  II),  also  indicate  rather  poor  fits  to  the  data,  as 

fv  2  2  2 

does  the  process,  ^  ^  [7]  (not  shown  on  Fig.  1.1). 
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While  the  tangent  plane  approach  of  [4],  as  used  here,  for  example, 
is  applicable  for  f  =  0(>  50  Hz),  it  is  not  suited  to  situations  where 
ko^c  not  ^ar9e’  c^-  (3. 1)- (3.3) .  We  must  then  use  the  Method  of 
Small  Perturbations  (MSP)  for  lower  frequencies  0(<  30  Hz)  where  the 
Rayleigh  parameter  is  small;  see  Secs.  9. 1-9. 6  of  [ 1 1 ] ,  esp.  Sec.  9.6 
vs.  Sec.  3  of  [5];  for  MSP  generally,  see  [12],  Chapter  3.  At  these 
low  frequencies,  the  coherent-scattered  energy,  in  the  specular  direc¬ 
tion,  dominates  the  incoherent  scattered  energy.  This  is  the  reverse 
of  the  "high-frequency"  or  large  Rayleigh  number  cases,  where  now  the 
surface  has  a  large  (vertical)  roughness  scale,  so  that  coherence  is 
destroyed,  even  in  the  specular  direction,  and  only  the  incoherent  scat¬ 
ter  is  observable,  cf.  (2.15 ) -[43  vs.  ( 2 . 13 ) - [4] .  In  our  present  study, 
as  noted  earlier,  we  shall  consider  only  the  intermediate  and  "high" 
frequency  situations,  when  kQ£  >(  >)  1  is  obeyed. 


Scattering  Cross  Sections:  Model  la 


Model  la,  cf.  Fig.  2.2b  and  Sec.  2.1,  consists  of  a  single  surface, 


which  we  arbitrarily  divide  into  low-  and  high-wave  number  portions, 
analogous  to  the  treatment  of  ocean  wave  surfaces  (Sec.  2.1,  p.  6,  and 
Eq.  2.1  of  [4]).  The  reason  for  this,  of  course,  is  that  the  low-wave 
number  part  provides  the  "facet"  or  "optical"  scatter  componiMs  while 
the  high-wave  number  portion  is  associated  with  the  Bragg  scatter  terms -- 
all  under  the  condition  here,  (3.3),  for  the  Kirchoff  or  tangent  plane 

method  [f  >  o(50  Hz)],  cf.  [13]. 

0  ~(n\ 

The  scattering  cross  sections  0:  '  for  the  incoherent  component 
of  reradiation  are  defined  and  discussed  in  Appendix  A. 2.1  (see  also 
Sec.  2.2  of  [4]).  For  Model  la,  (3.5),  we  have  here  the  first  two 
terms  of  (2. 13)- [4]  for  the  "high-frequency"  regimes,  i.e.,  large 
Rayleigh  numbers,  e.g.,  f  >0(1  kHz),  cf.  Table  3.1: 
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Our  general  conclusion  from  the  above  is  that  for  the  ice  profile  of 
the  Lyon  data  types  [1],  [5],  cf.  (2.9),  we  may  consider  the  tangent 
plane  or  Kirchoff  method  of  scattering  analysis  to  be  applicable  (with 
appropriate  corrections  0(1-3  db),  as  usual,  for  small  grazing  angles, 
cf.  remarks  pp.  3,  37,  44),  for  frequencies  above  0(50  Hz),  and  possibly 
above  30  Hz.  [For  convenience,  we  shall  designate  the  "intermediate 
frequency  range"  as  0(0.10  kHz  -  1  kHz)  here,  with  the  "high-frequency 
region"  defined  as  0{>  \  kHz).] 

The  Rayleigh  numbers,  R^,  for  these  surfaces  and  frequencies,  are 
typically  (for  a2  =  5.68  =  (2.38)2  (m2),  cf.  (2.6))  a -id  small  grazing 
angles,  here  0oT  =  80c  (for  which  cos  0qT  =  0.174): 

\  =  °-boko  =  (ko0;2  cos  6oT)2,  (backscatter)>  cf.  Eq.  (A. 2-12),  (3.4a) 
=  (0.828  kQ)2  ,  0oT  =  80°  ,  (3.4b) 

which  for  a  range  of  intermediate  frequencies,  0(30  Hz  -  1  kHz),  is 
shown  in  Table  3.1  below.  Frequencies  0(30-50  Hz)  are  borderline  for 
the  applicability  of  the  Kirchoff  method,  at  any  6oT  (<  90°).  We  note 
that  for  small  grazing  angles  (0^5*80°),  fQ  =  (30-100+  Hz)  are  "low- 
frequency"  from  the  viewpoint  of  the  Kirchoff  approach,  i.e.,  have  small 
Rayleigh  numbers,  while  f  =0(5*1  kHz)  are  "high-frequency,"  with  large 
Rayleigh  numbers.  At  larger  grazing  angles  (0oT<80°)  frequencies 
0(5*50  Hz)  are  all  "high-frequency"  in  this  sense. 


Table  3.1  Rayleigh  Numbers  for  Backscatterinq;  9  j  =  80 0 


fQ  (Hz):  (sec-1) 

ko  =  2tt/xo  ^ 

30 

50 

0.126 

1.09  •  10"2 

50 

30 

0.209 

3.00  •  10"2 

100 

15 

0.419 

1.20  •  10"1 

200 

7.5 

0.838 

4.82  •  10-1 

o 

o 

*3- 

3.75 

1.68 

1.93 

1  kHz 

1.5 

4.19 

12.0 

2  kHz 

0.75 

—  . . 

8.38 

48.1 
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(i).  The  Kirchoff  approach  requires  that  kQPc >  1;  or 

2ttpcA0  >(>)  1,  i.e.,  the  local  radius  of  surface  curvature 
must  be  large  compared  to  the  incident  wavelength,  XQ. 

This  is  equivalent  to  the  condition 

kQ£c  >(>)  1  ;  J l  =  (horizontal)  correlation  distance.  (3.1) 

(ii).  For  the  (isotropic)  bottom  ice  profile,  (2.10b)  gives  the 

normalized  directional  correlation  function  p  (Ar).  From 

tables  of  Ki(x),  p.  737,  ClOl,  it  is  easily  found  that 

xK^(x)  =  e"1  at  x  =  x_^  =  1.65  =  k(;Ar_1,  so  that,  for 

Ar  ,  =  £  we  have 
-1  c 

1.65(kQ/kc)  >(>)  1,  or  (Xc/Ao)(1.65)  >(>)  1  (3.2) 

as  a  defining  condition  that  the  Kirchoff  method  is  valid. 
Since,  empirically,  kc  =  0.05,  we  require 

H>)  x0>  2.07  ‘  10Z  >(>)  *Q.  (3.3) 

(ill).  For  f  =  30  Hz,.*.  Xq  =  50  m,  and.*.  2.07  •  102  >(>)  50, 

which  is  not  a  very  strong  inequality.  However,  for  fQ  = 

50  Hz,  Xq  =  30  m  and  2.07  •  102  >(>)  30,  which  is  better. 

At  fQ  =  100  Hz,  Xq  =  15  m,  and  the  inequality  is  -14:1, 

which  is  sufficient.  Moreover,  defining  i  according  to 
_  1  ^ 
the  e  -scale,  cf.  (3.2),  is  quite  strict:  after  PJl  )  = 

- 1  c 

e  =  0.368  there  is  still  noticeable  (radial)  correlation 

left.  We  may  therefore. choose  a  distance  such  that 
st<Voi>  * 0  .01,  i.e.,  essentially  negligible  surface  cor¬ 
relation,  for  which  k^  ±  6.5  (on  using  [10]),  so  that 
(3.3)  is  modified  to  2ir(6 . 5)/Q .05  =  8-16  •  102  >(>)  \  ,  which 
in  turn  permits  a  lower  frequency,  say  20  Hz  (or  XQ  =  75  m), 
for  8.16  •  102/75  *  11  >(>)  1. 
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where  again  (2.14a,b)  apply.  Now,  however,  does  not  enter  explicitly, 
cf.  (2.15).  The  gaussian  processes  described  statistically  by 

(2.16a,b),  still  apply,  as  do  the  conditions  (2.17).  The  statistics  of 
noted  in  (2.18)  must  once  more  be  evaluated,  in  any  case. 


2.4  A  Model  for  the  Marginal  Ice  Zone  (Mil)  (Model  lib) 

Model  lib  above  is  also  a  reasonable  candidate  to  describe  the  com¬ 
bined  presence  of  ice  and  open  water  which  is  typical  of  the  Marginal 
Ice  Zone  (MIZ).  Now  we  replace  the  keels  (in  AK  above  (2.19))  by  the 
(plate)  ice,  e.g.,  r  e  Aj,  and  5$  *  Cocea„  Surface,s. 

£  e  A- A j  t  Aj,  e.g. , 


C(r. t)  =  -^2CI(r,t) • 

=  ~—z  Hte'V  * 


£  e  Aj;  1 

r  £  A-Aj  t  Aj  )  , 


(2.20) 


where  Cj  =  Cj(r)  if  there  is  no  relative  motion  of  the  ice  vis-a-vis 
the  platform  in  question  (cf.  (2.8)  et  seq.).  Now  £<.(£, t )  represents 
the  (moving)  ocean  wave  surface,  as  before,  which  is  subject  to  wave 
generation  by  the  wind,  in  the  usual  way  [vide  Sec.  1,  Sec.  2.1,  of 
[4]]. 


The  further  development  of  the  analysis  for  Models  IIa,b,  including 
the  desired  scattering  cross  sections,  is  reserved  to  a  subsequent  re- 

e  with  Model  I  only. 


3.  Scattering  Cross  Sections  for  Model  la 

Here  we  make  direct  use  of  the  results  and  derivations  of  the  au¬ 
thor's  recent  study  [4].  We  shall  accordingly  refer  to  specific  rela¬ 
tions  therein  by  equation  number,  with  -[  ]  suffixed,  e.g.,  (2.13)  of 
[4]  becomes  (2.13)-[4],  etc. 

Since  our  analysis  is  based  on  the  Kirchoff  or  tangent  plane  ap¬ 
proach,  the  results  are  mainly  applicable  in  the  intermediate  through 
"high  frequency"  ranges,  which  may  be  numerically  established  here  from 
the  following  argument: 
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posed  processes,  then  W^(k)  is  a  sum  of  two  separate  processes,  so  that 
and  p^  may  overlap  (in  wave-number)  space.  This  partition  is  then 
made  on  the  basis  of  the  relative  contributions  to  the  total  process, 
i.e.,  on  an  intensity  basis,  as  indicated  empirically  by  the  parameters 
as,  aK  in  (2.2)  and  (2.16a,b): 

(2) .  cK,  cs  are  independent,  by  the  argument  of  Sec.  1  above;  j 

(3) .  <5k  is  independent  of  (and**  (5x,y^K:  the  ^ of 

the  keel  is  independent  of  keel  elevation,  as  well  as  of  the! 

small er-scale  surface  component,  . 

(4) .  6k  can  be:  (i),  quasi -periodically  distributed  [in  the  ( 

(x1 ,y' )-plane];  or  1 

(ii),  entirely  randomly  distributed  in  this  plane,\ 
so  that  overlapping  keels  are  possible.  J 
In  the  analysis  leading  to  scattering  cross-sections  (and  spectra  where 
appropriate,  i.e.,  moving  platforms),  we  shall  need  the  following  sta¬ 
tistics  for  <5^: 

<V  i  <SK1SK2>  * 

again  with  the  usual  assumptions  of  homogeneity  and  isotropy,  cf.  (2.7) 
et  seq.  From  (2.14b)  it  follows  that  m$|max  =  1  =  m6(0,0),  at  least, 
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where  A^  is  the  region  occupied  by  the  keels,  cf.  Fig.  (2.2a).  [The 
positive  z-direction  is  upward  from  the  reference  (x,y)-  and  parallel 
(x1  ,y‘  )-planes ,  with  £  =  j^x  +  J^y  =  j_xx'  +  j^y' ,  etc.  ]  In  terms  of 
the  subsequently  required  integrals  over  the  domains  of  acoustic  "illu¬ 
mination"  of  the  underice  surface  we  have 


jr  e  Ak  J  $K(r,t)(  )dA  ;  r  e  (A-A^)  -*■  /^[l  -  6K(X»t)  ](  )dA,  (2.14a) 


with 


V1’  rtAK; 


-  0,  r  £  A 

—  I\ 


(2.14b) 


Here  n„  is  the  unit  normal  to  the  ice-keel  surface  component,  given  by 


nK  =  !(ixcx  +  Vy  '  ^ +  5x  +  cv  tic  ; 


y  »K 


34 


(2.15) 


=  -j2  when  r  e  A^  4X  =  ^  ,  etc. 


The  following  assumptions  and  conditions  are  next  imposed: 


(1).  4k,  4s,  and-  4 Kx ,  4^,  are  gauss  processes,  now  with  the  c.f.'s 

-im$5-a^2/2  -ic^ -ic2m$  -  (o2/2)[Cj+C2+2C152Ps] 

(2.16a) 

•  m  r  „2r2/9  -iC1m^-U9m„-(o2/2)[c?+C9+2C15,P^] 

FjdOjc  =  e-imK5’°KC  /2;  F2(44r*6^^“-J-X— ^ - L-L-l  2  K 


(2. 16b, 


where  m^,  a^,  mK,  are  given  (in  the  case  of  Lyon's  data)  by  (2.3), 
and  we  have  taken  account  of  the  fact  that  the  underice  elevations  are 
negative  in  our  reference  system,  cf.  Fig.  (2.2a).  Here  p^  and  pk  are 
normalized  directional  covariances  of  these  elevations.  From  the  iso¬ 
tropic  assumptions  of  the  empirical  fit  to  the  Lyon's  profile  given  by 
Mellen  et  al.  [5],  and  p^  are  to  be  obtained  by  some  suitable  parti¬ 
tion  of  Wj(k) ,  (2.9),  or  W2(k,4>),  (2.11b).  If  we  follow  the  underlying 
assumption  of  Model  I  la  above  here  that  the  S-  and  K-surfaces  are  super- 
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Model  1 1  b 


Model  I  la 


Figure  2.2a:  Model  II 


Model  II :  Two-dimensional  underice  profiles  (c),  Eq . 
(2.13) ,  with  distinct  keels  (CK);  K  are  the  (random) 
regions,  projected  onto  the  (x,y)-  or  (x',y':ca0)  planes 
by  the  specific  keel  structure.  Model  I la:  Cc  rides  on 
Model  lib:  and  are  spatially  separate. 
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(gaussian)  underice  profiles,  C,  of  Model  I,  where  keel  structure,  as 
such,  is  implicitly  subsumed  in  the  profile  statistics,  cf.  Fig.  1.1. 

Two  variants  of  Model  I  will  be  examined:  Model  la,  which  is  a  single¬ 
surface  model,  is  the  simplest  choice,  and  one  which  we  shall  find  needs 
to  be  extended,  to  account  for  observed  scattering  cross  sections.  This 
extension  is  Model  lb,  which  now  includes  an  additional  small-scale  rough¬ 
ness  component,  riding  on  the  single-surface. 

Before  we  continue  on  to  the  results  for  scattering  cross  sections 
in  Sec.  3  ff.  for  Model  I,  let  us  concisely  describe  Model  II,  which  is 
a  viable  alternative  to  Model  I,  at  least  in  those  interactions  where  a 
keel  structure  can  be  identified.  As  in  the  case  of  Model  I,  two  vari¬ 
ants  of  Model  II  appear  plausible:  Model  Ila,  which  is  characterized 
by  a  comparatively  small-scale  (small  m,  a)  elevation  process  which 
rides  upon  the  (quasi-)  random  keel  structure  (large  in);  and  Model  lib, 
where  the  "small-scale"  profile  and  the  keel  profiles  occur  separately 
in  space.  This  latter  form  also  appears  reasonable  for  modeling  the 
Marginal  Ice  Zone  (MIZ)  cases. 

2.2  Model  Ila:  A  Two-Component  (Two-Dimensional 


Unden ce  Profile,  with  Keels 


For  this  model  we  postulate  that  there  is  a  distinct  keel  structure, 
with  profile  component  upon  which  "rides"  the  smaller-scale  bottom 
irregularities,  whose  profile  component  is  1n  the  manner  of  Fig. 

2.2a.  This  is  to  be  contrasted  with  our  unstructured  two-component  model 
(Sec.  2.1),  where  no  explicit  keel  structure  is  distinguished  but  is  re¬ 
garded  as  part  of  a  single  profile,  C,  governed  by  (2.1)  etc.,  as  sketched 
in  Fig.  2.2b. 

Accordingly,  the  underice  profile  for  our  Model  II  can  be  expressed 


^(r.t)  =  -iz^K(r,t)  +  nK(r,t)^s(r,t) ,  I  e  Ak 


(2.13) 


*  -Vst'1’ 


*The  time-dependence  is  indicated  explicitly,  to  include  possible 
platform  motion,  cf.  (2.8)  et  seq. 
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from  Appendix  I,  where  K^(x)  is  a  modificed  Bessel  function  of  the  second 
kind.  The  corresponding  transform  relation  is 

oo 

W,(k)  =  a2  J  eikAr  p  (Ar)d(Ar)  =  Eq.  (2.9),  etc.  (2.10c) 

**■  ^  ^  _  CO  ** 

The  associated  two-dimensional  wave  number  intensity  spectrum  for  this 
isotropic  model  is  accordingly  given  by 

(2.11a) 
(2.11b) 


in  which 

00  00  - 

2  _  r"  d(j)  /  W_(k,4))kdk  =  2  /  W.(k)dk  =  C J\r  <  *  (2.11c) 

5  0  0  2  0  1  1  c 

2 

cf.  Eq<  (4)  of  [5],  and  Appendix  I.  We  note  that  while  a ^  can  be  ana¬ 
tomized  into  the  component  contributions  of  cf.  (2.5),  (2.6), 

it  is  not  possible  to  do  this  for  the  covariance  or  wave  number  spectrum 
(2.8)-(2.11)  in  our  present  Model  I. 

2 

For  these  empirical  models  (2.9),  (2.11)  clearly  a ^  is  finite,  as 
it  should  be.  However,  we  note  that  the  corresponding  wave  number  spec¬ 
tra  of  the  profile  slopes,  c  ,  c  ,  or  3e/3x,  etc.,  are 

a  y 


Wl(k)3?/3x  =  ’  tyM)^  =  k2W2(k,4>) 

x-y 


(2.12) 


2 

so  that  the  mean-square  slopes  a  are  logarithmical ly  divergent.  Of 

sX  yj 

course,  physically  this  is  not  acceptable,  and  we  must  therefore  require 
a  faster  fall-off  of  the  model  W^  ^  with  wave  number,  at  the  high  fre¬ 
quencies.  As  yet  there  do  not  appear  available  data  which  allow  us  to 
extend  the  results  of  Fig.  (2.1)  beyond  Xg <  6  m.  Consequently,  we  are 
unable  at  this  time  to  determine  a  ;  (see,  however,  Sec.  3  ff.). 

On  the  whole,  nevertheless,  with  (2.5),  (2.6),  (2.9),  (2.10)  in 
(2.4)  for  Fj(i£)^  and  (2.7)  for  F2^  we  are  specifically  equipped  to  use 
the  author's  recently  developed  methods  and  results  [4],  (1984)  formally 


to  obtain  the  desired  scattering  cross-sections  for  the  composite 
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The  new  quantity  p  is  the  normalized  covariance  of  the  (possibly) 
relatively  moving)  underice  surface  elevation,  defined  by 


pr(£r.T=!£rl/v  ) 


^ '£2 ’ t2 ^ 


{5(!l1»t1)C(r1-»Ar,t1+ t)  -  C(r,t)  ^)/o2; 


Ar  =  j;2-_r1;  x  =  tg-tj  , 


(2.8) 


where  £  =  lxx  +  j^y  is  a  vector  in  the  (x,y) -reference  planes,  cf. 

Fig.  (2.2a)  ff.,  and  vp  is  the  (uniform)  speed  of  the  reference  plat¬ 
form.  Thus,  if  there  is  platform  motion,  the  (random)  ice  profile 
appears  to  be  moving,  like  a  random  wave  surface.*  (With  no  relative 
motion  p„(Ar,r) -> p_(Ar,0)  =  p_(Ar),  as  expected.)  Mote,  as  required, 
when  |Ar|-+«,  p^  +  0:  there  is  no  "dc" -component  of  ?  in  p?,  and  accordingly 
F2 ! p-H3~Fl^i^l)Fl(i^2)’  where  Fi  is  not  given  anymore  by  (2.4);  cf.  (2.6)  et  seq 
Next,  we  use  Mellen  et  al.'s  [5]  isotropic  empirical  fit  to  the 
one-dimensional  wave  number  (intensity)  spectrum  of  the  profile,  Eq.  (2), 
[5],  obtained  by  FFT  analysis  of  the  profile  shown  in  Fig.  1.1  with 
"dc"  (or  mean)  component  removed.  This  spectrum  is  (See  Appendix  I  ff.): 


I 

j.W1(k)5  =  trC1/(k^  +  k2)3/2;  C1  =  cr2k2  .  0<013;  ^  =  0.05  nf1 


(or  Xc  =  126  m). 


(2.9) 


which 


is  given  in  Figure  2.1.  The  standard  deviation  a  ,  of  found 


from  the  analysis  of  these  data  is  -  2.3  m,  cf.  remarks  after  Eq. 

(2)  of  [5],  a  result  which  is  in  excellent  agreement  with  that  of  our 
composite  model  above,  (2.6),  namely,  a =  2.38  =  2.4-  (m).  The  asso¬ 
ciated  (normalized)  covariance  function  for  (2.9)  is  (vp  =  0): 

00  . 

P?(£>  =  h  /  dkwi(k)^"ik|^r|/V»  I  Ar  I  ■  lr2-ri !  •  Ar,  (2.10 


tt  L  Wl(k)  cos  k  Ar  dk/a^  =  kcK1(kcAr)Ar, 


(2.10 


*For  random  platform  movement,  in  the  general  "bi-static"  case  of 
separate  moving  source  and  receiver  platforms ,  see  the  author's  1977 
paper,  [9]. 
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From  well-known  results  ([6],  Eqs.  (1.30),  (1.31))  we  have  the 
associated  characteristic  function  (c.f.) 


F 1  ( i  £ )  =  /  w1(?)eK^dC  =  t^e 


ImjC-ofC V2.  imKC-a^C  /2 

+  b2e 


with  the  following  first  and  second  moments 

1  =  +  b2mK  =  5.1  (m);  C2  =  b^a2  +  m2)  +  b2(a2  +  m2)  =  31.3  (m2) 

(2.5) 


2  -2 

=  ?  -C  =  var  s 

=  b^c2  +  b2<32  +  b^m2  +  b2m2  -  (b^  +  b2mK)2 

2  2  2 
=  b^  +  b2a^  +  bib2^rnS“mK^ 

=  5.68  (m2)  =  ( 2 . 38 ) 2  (m2) 


(2.6) 


from  (2.2).  We  note  that  if  we  assumed  that  the  underice  elevation  c 
were  a  composite  (gaussian)  sum,  e.g. ,  c  =  bjC$  +  b2^«>  then  is  cor~ 
rectly  given  in  (2.5),  but  then  £2  f  ^^omp>  cf.  (A. 3-13).  In  fact,  we 
show  in  Appendix  A. 3  that  (I),  it  is  not  possible  to  construct  a  second- 
order  gauss-composite  process  which  reduces  to  (2.1)  or  (2.4). 

On  the  other  hand,  we  need  primarily  the  second-order  pdf  (or  c.f.) 
of  the  profile  ?,  in  order  to  determine  the  desired  (incoherent)  back- 
scatter  cross  section.  In  view  of  (I),  Appendix  A. 3,  remarked  on  above,  we 
must  restructure  our  gauss-composite  model  to  include  one  which  has  a  (still 
gaussian)  second-order  pdf  (and  c.f.),  again  because  of  the  Central  Limit 
theorem  argument  cited  in  Sec.  1.  Accordingly,  we  choose  the  composite 
gauss  process  s  =  b^  +  ^K’  and  requ^re  ^  t0  have  the  same  first-order 
first  and  second  moments  as  (2.1)  above,  cf.  (2.5),  (2.6),  accepting  the 
rather  poor  first-order  fit  to  the  data,  cf.  Fig.  2.1.  [For  the  coherent 
scatter  cross  sections  we  would  use  the  well-fitted  results  (2.1)  and  (2.4) 
directly.]  Consequently,  the  following  second-order  characteristic  func¬ 
tion,  from  Eq.  (7.10)  of  [6],  can  be  written 


^1  ’^2^c^  ”  ® 


i51C+iC24-c2[C2+?2+2E1C2Pc]/2 


U  =  blCs+b2cK), 


where  we  have  postulated  homogeneity  (and  stationarity)  for  our  present 
_  2 

model.  Here  ?  and  a_  are  given  by  (2.5),  (2.6). 
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2._  Underice  Surface  Elevation  Models I  and  II 

Here  we  lay  the  statistical  foundation  for  both  our  Models  I  and  II, 
referred  to  above  in  Sec.  1,  although  in  this  Report  we  shall  use  only 
Model  I  in  carrying  out  our  determination  of  the  desired  scattering  cross 
sections.  Accordingly,  let  us  begin  with  Model  I. 

2.1  Model  I:  An  Unstructured  Two-Component  Underice  Profile 


From  our  fit  to  the  Lyon  histogram  illustrated  in  Figure  .1.1  we 
obtain  the  following  first-order  pdf*  for  the  underice  "profile"  or  sur 
face  elevation  : 


U-mc)  /2a 


(C-mJ  /2a 


where  we  find  that,  specifically, 
a. 

b^  *  TTT-  =  0.93 ;  b,  =  *  0.070;  a.  =  0.22;  =  0.016; 

l  a.+a,,  2  a.+a«  —  1  2 

1  c  L  c  and  b]+bz  =  1.00; 


mS  =  4.50  m 


mK  =  12.5  m 


Oj  *  1.00  ml  ’  a^  =  3.0  m 


It  is  clear  that 


/  w.(C)dC  *  /  w.(C)dC  =  1 
0  1  -®  1 


here,  in  view  of  the  large  means  (m^,  m^)  and  comparatively  small  vari¬ 
ances,  so  that  the  "tails"  of  both  the  fitted  gaussian  densities  are 
quite  negligible  for  c<0.  Thus,  the  underice  "profile,"  or  elevation 
C,  is  described  by  the  weighted  sum  of  two  gaussian  distribution  (den¬ 


sities)  (which  is,  of  course,  not  the  same  as  the  composite  gauss  process 
consisting  of  the  weighted  sums  of  the  separate  elevations  zs>  5^, 
cf.  Appendix  A. 3.)  We  call  z  here,  cf.  (2.1),  a  gauss-composite  process. 


♦For  the  moment,  in  Model  I,  we  regard  z  as  positive  downward;  since 
z  is  positive  upward,  cf.  Fig.  (2.2a),  in  the  platform-underice  geometry 
of  Section  3  ff . ,  we  will  accordingly  change  the  signs  of  5,  etc.  appropriately 
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Because  of  the  large  Rayleigh  numbers  involved,  theories  based  solely 
on  Bragg  scatter  are  inadequate,  except  at  very  low  frequencies  and/or 
large  grazing  angles.  Consequently,  we  must  develop  the  first-  and 
second-order  characteristic  functions  of  the  profile  elevation  C.  Then, 
we  must  artificially  decompose  the  single  elevation  process  C  (Model  la) 
into  two  further  ( independent)  processes  at  some  sufficiently  short  sur¬ 
face  wavelength,  or  corresponding  wave  number  (=2Tr/Xp),  so  as  to  ob¬ 

tain  the  scattering  contributions  of  the  "rough1  (i  .e. ,  Kirchoff  "facets") 
and  the  "smooth"  (or  Bragg)  surface  components,  just  as  is  done  in  the 
conventional  treatment  of  wind-wave  surfaces  as  distinct  gravity  plus 
capillary  wave  surfaces  [8].  (A  similar  approach  is  used  for  the  large- 
scale  surface  component  of  Model  lb  ff.) 

The  specific  dichotomy  depends,  of  course,  on  the  wave  number  (in¬ 
tensity)  spectrum  of  the  underice  surface  (cf.  Eq.  2.3  and  Sec.  7.3C  of 
[4]),  examples  of  which  are  modeled  in  [2],  [5]  (Sec.  2)  from  empirical 
data.  Thus,  our  theoretical  approach  in  this  work  is  a  direct  extension 
and  modification  of  our  earlier  results  [4],  which  provide  the  substan¬ 
tive  analytical  foundation  of  the  present  effort.  We  shall  for  the  most 
part  be  content  to  present  the  needed  final  results,  referring  the  reader 
to  C4]  for  the  detailed  background.  What  is  new  here  vis-a-vis  earlier 
work  [for  example,  [2],  [5],  and  refs,  therein]  is  (1),  explicit  statis¬ 
tical-physical  models  of  the  underice  profile,  which  permit  an  incoherent 
as  well  as  coherent  scattering  analysis  at  intermediate  and  high  fre¬ 
quencies;  and  (2),  the  specific  results  of  that  analysis,  albeit  limited 
to  Model  I  here.  Model  I  itself  is  simplified  basically  to  a  single 
gaussian  process,  5,  now  with  means  and  variances  given  by  those  of  the 
first-order  gauss-composite  fit  to  the  Lyon  data  [cf.  (2.1),  (2.2)  et  seq.] 

Finally,  this  Report  is  organized  as  follows:  Section  2  gives  an 
analytic  description  of  the  aforementioned  two  profile  models  (I,  II), 
with  their  pertinent  statistics.  Section  3  provides  backscattering 
cross  sections  of  Model  la,  and  Section  4  includes  the  extension  of 
Model  la  to  Model  lb,  to  provide  a  possible  mechanism  whereby  agreement 
between  theory  and  experimental  data  can  be  effected.  (A  full  treatment 
is  reserved  to  a  subsequent  report.)  Section  5  concludes  our  analysis 
with  a  concise  summary  of  the  principal  results  and  conditions  under 
which  they  apply.  [Appendixes  A.1-A.3  provide  appropriate  special 
results  and  analyses.] 


4 
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.  J_  p2a of  t  <2°oy)21 

8(0)  ^2  i  /‘‘a«*‘‘oy*°oz  \  2  e  2bo  _  j>  °Ly 

i„c  Rt»l  -  So  *Z  ^  o02/2  } 


+  X  N<°> 


(3.6) 


^■Ni  w  (2a  )  •  Wu(2a.  .  k  )  [ 

LH  i  nc  H  ox  o7  j  ^ 


The  terms  in  {  }  represent  the  results  of  the  "physical -optics"  cum 
"perturbational "  approach,  i.e.,  here  respectively  the  facet  and  Bragg 
scatter  terms,  cf.  Sec.  7.3  of  [4].  The  former  is  seen  to  be  independent 

of  frequency,  while  the  latter  is  a  first-order  Bragg  component,  since 

2 

a  <<  1,  or  more  precisely,  RUr«  1,  with  small  correlations.* 

H  /M  “4 

In  (3.6),  N'  ’  "  +  -i  1  <---P  11  ''"KnJwinn  -M,  = +• 


LH-inc 


is  a  "tilt-factor,"  embodying  the  fact  that  the 


high-wave  number  surface  component,  cu,  is  tilted,  i.e.,  presents  on 
the  average  a  larger  grazing  angle  4>(=  tt/2  -  6qT)  toward  the  plane  of 
incident  radiation.  The  quantity  WH  is  the  wave  numner  intensity  spec¬ 
trum  of  into  which  the  (single)  profile  4,  (3.5),  is  dichotomized, 

on  2 

while  aLx,  aLy  are  the  variances  of  the  (x,y)-slopes  of  the  ”low,,-wave 
number  component,  Cp,  of  the  underice  surface  in  this  model,  cf.  (3.5). 
(The  symbol  1  denotes  (x,y)-plane  or  (x,y) -components  only,  and  a  , 

U  A 

etc.,  are  differences  of  the  components  of  the  incident  and  reflected 
waves'  unit  vectors,  -  1r>  ^or  the  basic  one-reflection  geometry  of 
this  study,  viz.  Figs.  A. 2-1,  A. 2-2;  bQ  =  cos  9qT  +  cos  0oR.  For  a  de¬ 
tailed  description  of  the  various  elements  of  (3.6),  and  (3.7),  (3.8) 

following,  see  Appendix  A. 2,  including  the  reflection  and  shadowing 
?  ~ 5" 

factors  (i?g,  S  )  in  Sec.  3.2  below. 

For  the  coherent  cases  the  scattering  cross  section  is  defined  by 
(A. 2-5)  and  is  found  for  all  Rayleigh  numbers,  R;-(>  0)  such  that  the 
Kirchoff  condition  (3.1)- (3. 3)  is  obeyed,  here  for  f  (>  o(50  Hz)),  to  be 


*This  assumes  a  "splitting"  wave  number  kp  large  enough  that 
00 

°H  WH(J^/(27T)2  <<  1  (m2)>  cf.  Sec.  3.3C,  Eq.  3.21)  ff. 
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d(0)  _ 

coh 


i  it 

O  0 

- r 

16ti 


<°lkl  -*yox/l  *  a 

e  •  e 


-*•  0,  "high-frequen".  ;es ,"  R  »1, 


here  for  f  =  0{>  1  kHz),  cf.  Table  3.1  above:  very  rough  surfaces, 
e.g.,  R^  »1,  destroy  coherence,  as  expected.  Lin  addition,  since  the 
beam  parameters  A,B  «1,  cf.  (6.32)-[4],  the  second  exponential  also 
ensures  the  vanishing  of  3cQh,  unless  the  Snell  angle  (aQx  =  a  =  o) 
is  chosen.]  Note  that  this  coherent-scatter  cross  section  depends  on 
the  area  illuminated  (via  4^),  as  distinct  from  the  incoherent  cases, 
(3.6),  which  are  always  area-independent,  subject,  of  course,  to  the 
conditions  cited  in  (A. 2-15).  Our  results  (3.6),  (3.7)  hold  for  ar¬ 
bitrary  directions  of  illumination  and  observation. 

Various  particular  cases  of  (3.6),  (3.7)  can  be  obtained  directly 
by  appropriate  specialization,  e.g.,  "high"-  and  "low "-frequency  back- 
and  forward-(specular)  scatter,  etc.,  from  the  author's  results,  pp. 
18-24  of  [4]  (always  subject  to  the  Kirchoff  condition  (3.3)-(3.4) ) . 
Since  we  shall  show  presently  [Sec.  3.3]  that  Model  la,  (3.5),  is  ap¬ 
parently  incapable  of  acceptable  agreement  with  empirically  measured 
cross  sections,  we  shall  limit  our  treatment  to  the  special  cases  of 
backscatter,  for  which  (3.6)  reduces,  on  setting  4>qT  *  iT/2,  cf.  Fig. 

A-2'1,  6oR  =  9oT’  t0 


e(°) 

inc 


back 


*  «0  52 


Rr»l 

6^-/2 


[ 


-  tan2eoT/2^y 


8lT0Lx°Ly  cos  9oT 


16tt 


°INLS)(6oT)1nc“H<0'2koslr 


A 


t 


e(o) 

°inc-L 


,(0) 


inc-H  )R^»1 


(3.8) 

(3.8a 


Again,  the  first  term  of  (3.8)  is  independent  of  frequency,  in  this 
"high-frequency"  (R^  »1)  approximation,  which  as  noted  above  is  the 
"specular-point,"  facet,  or  geometrical  acoustics  solutional  form  of 
the  "full-wave"  approach  of  Bahar  and  Barrick  et  al.  [13],  [14],  [15], 
(and  also  references  [5],  [22],  [24],  [31],  [44]  in  [4]).  The  second 


term  of  (3.8)  embodies  the  (first-order)  Bragg  scatter  associated  with 

the  small-scale  surface  component  qu.  Whereas  the  specular-point  con- 
(0)  M 

tribution,  is  the  dominant  part  of  the  cross  section  for  moder¬ 

ate  to  small  9 0y,  it  vanishes  rapidly  for  small  grazing  angles  (9qT 
large),  so  that  in  this  surface  model  only  the  perturbational  Bragg 
term,  3  u,  contributes.  We  shall  use  this  observation  in  Section 
3.3  below  to  test  the  applicability  of  our  Model  la,  (3.5). 


3.2  The  Reflection,  Tilt-,  and  Shadowing  Factors: 

~p2  m  ( ® )  c2  . 

V  TH-inc’  5  • 


With  wind-wave  surfaces,  which  are  pressure  release  surfaces, 

r/(pc0)air  ~  08 •  As  a  result  the  (local  plane-wave)  reflection 
coefficient  R  i s  -  -1,  cf.  (3.9a)  ff.  Moreover,  since  the  rms  wave 


(pco^wate 


slopes  are  small  0(<  5*  10"2),  the  tilt-factor  (N^/l£>  =  N)^)  is  li 


( 0 )  /  1  c  =  m  ( 0 ) 


ike- 


wise  small  °U0~2),  cf.  (3.2d,e)-[4],  and  unless  0oj>  850,  the  shadowing 


coefficient  s1*  is  always  close  to  unity.  With  ice  surfaces,  however, 
these  perturbational  conditions  do  not  fully  hold:  as  we  shall  see  from  (3.9) 


below,  at  these  frequencies  o£*50  Hz),  RQ  is  +1,  while  quite  large  slopes 
a2(x=y)=o(<1.5)  may  occur.  This,  in  turn,  generates  larger  tilt-factors,  as 


we  would  expect,  and  at  the  same  time  produces  more  shadowing.  Consequently, 
we  must  reexamine  our  earlier  evaluations  of  both  the  tilt-factors  and 
shadowing  coefficients  used  in  our  wave  surface  analyses  [4],  to  take 
these  modifications  into  account. 

Accordingly,  let  us  begin  first  with  the  reflection  coefficient: 


A.  The  (Plane-Wave)  Reflection  Coefficient,  R 

For  our  far-field  application  here  the  local  plane-wave  reflection 
coefficient  is  a  reasonable  approximation,  even  when  averaged  and 
removed  from  under  the  integral  sign  in  the  Total  Surface  Spreading 
Function  (TSSF),  cf.  (A. 2. 7a),  as  in  (3.6)-(3.8).  Specifically,  we 
have 


PTCTC0S  9  T  ■  Pu,cu,c°s  0«T 
/?  f 9  \  —  i  i  oT  w  w  or 

o'  or  '  PjCjCOS  e„T  +  Pj^cos  9' 


JoJ 


w  w 


oT 


where  (pc)T  -  2  for  ice  and  (pc) 


1  for  water,  and  9'T  =  angle  of 
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transmission  in  the  ice.  Since  at  these  frequencies  (perfect)  internal- 
reflection  occurs  for  ice,  Q^y  =  tt/2,  and 


2cos  0oj-O 
Ho  "  2cos  9oT+0  "  1 


(3.9a; 


[For  air-water  interfaces  { pc: ) j^.a r <<  1 »  and  so  RQ  *  -1,  as  usual.]  Thus, 
7?2  is  our  scattering  cross  section  here,  and  is  given  by 


rc  ±  1  =  0  db  . 
0  ice 


(3.10) 


B.  “Tilt-Factors"  N<°l1ne,  S<°>coh 

From  (8.28a)-[4],  and  from  (7.22)-[4],  (7.24)-[4],  respectively, 
for  the  incoherent  and  coherent  cases,  we  can  write  the  "tilt  factors" 
for  Model  la  here  as 

NLH-inc  =  ^2aoz^  ^2aox5x  +  2aoy?y  "  2aoz^  /£  *  ;x  =  3x*  etc‘  (3‘11^ 

For  the  symmetrically  distributed  slopes  assumed  here  (about  zT  =  if?  =  0), 

X  y 

(3.11)  reduces  directly  to 


NLH-inc  =  (2aoz)2  {(2aox)2aLx  +  {2aoy)Z°iy  +  (2aoz)2  } 


Lx,y  '  gx,y  ’ 


Choosing  <p  T  -  tt/2,  cf.  Figs.  A. 2-1, 2,  we  get,  cf.  (A. 2-11), 


(3.12a 


NLH-inc  ^oT=tt/2  ~  ^2ctoz^  ^2aoy^  aLy  +  ^2ao 3* 


(3.12b 


Similarly,  from  ( 7 . 22 )- [4] ,  etc.  we  have 


m(0)  _/^2aox*»x  +  2aoyfy  '  2ctoz^ 

LH-coh  " \  ,  2  2 

X  1  +  C  +  C 


(6c3a  ^+6t,^a  ^+a  ^) 
n  /  X  ox  y  oy  oz  \ 

8aoz  \  2  2/ 

\  1  +  C  +  C4  x 


taking  advantage  of  the  assumed  symmetrical  distribution  of  the  slopes 
(where  <4^>  =  <4  >  =  0,  etc.). 

X  X  r\ 

Inasmuch  as  <4>^  >  =  mean-square  slope  angles  of  the  profile,  we 
x  ,y 

have  at  once  by  definition 


Unlike  surface  waves,  where  =  (5.42*  10  t)  or  <4>t>"2  =  13.1°, 

lx  y 

cf.  (3.2d)-[4J,  for  20  knot  (mean)  wind  speeds,  we  may  expect  the  slopes 
of  the  bottom  ice  profile  to  be  considerably  larger.  This  is  borne  out 

by  (albeit  limited)  empirical  data,  which  suggests  that  > 2  =  0(30°), 

2-1  2  2 

or  °i  v  =  3.33  *  10  ,  [16].  Consequently,  we  cannot  replace  1  +  4+4 

LX,y  x  y 

by  unity  in  (3.13),  but  must  evaluate  the  averages  exactly.  We  continue 

to  assume  isotropic  bottom  surfaces  [cf.  Sec.  2],  e.g.,  =  o 

Accordingly,  since  4  is  postulated  to  be  gaussian,  in  accord  with  our 

Model  I,  cf.  Sec.  2  above,  so  also  are  4V,  4.,,  which  are  likewise  inde- 

x  y 

pendent  of  each  other  and  of  4.  In  fact,  we  can  write  explicitly 

-  ?iWy  ~^Lx+?Ly^2oL(x=y) 

W1^lv'^lu)  =  „  =  9 


r  Lx  Ly‘ 


2™LxaLy 


5L(x=y) 


2  2 

with  aLx  ^  =  cLx  .  Using  (3.15)  in  (3.13)  allows  us  to  write 
^LH-coh  5  -Kz  t(6olox  *  6a2y)I<2V)«.2z>  • 


isotropic 

(3.1 


where  (cf.  pp.  107,  108,  [4]) 

_  2  2  -(x2+y2)/2o2  .  2 

.(2)<a2M/  J  i£  or.Ae  -  ^dy, 

1  +  x^  +  y  2tto  0  1  y 

a2  =  1/2 a2. 


=  i  [1  -  a2  /  dy>  *  J  (1  +  a2ea  Ei(-a2))  , 
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since 


00  -a  z 
■  e 


r  =_ 
o  1  + 


;a  Ei  (-a2) ;  with  -ea  Ei  ( -a2 )  =  a-2  l  (-l)m  m!a'2m,  a2»l, 

m=0 

(3.17a) 


and  -Ei(-x)  =  Ei(x)  is  tabulated  in  [17],  [18].  Table  3.2  and  Fig.  3-1 

lQ\  ~(0) 

below  show  N^H-inc  and  f\H-coh  for  a  ran9e  ms  P^of i  1  e  slopes. 

2 

C.  Shadowing  Factor,  S 

We  shall  limit  ourselves  here  to  the  incoherent  backscatter  cases 
(ROT)  only,  for  this  study.  [A  summary  of  pertinent  results  for  and 
S  (ROT)  is  given  in  C_,  Sec.  7.4,  [4],  cf.  p.  109.]  Here  ^ack  =  Q(^  =  4>), 
Fig.  59  of  [12],  where  <J>(  =  tt/2  -  6qT)  is  the  grazing  angle,  as  before. 

The  pertinent  parameter  here  is 


a  -  (cot  0oTVaL(x=y)  (  =  a  of  Fi9-  59,  [12]),  \ 

■(tan^°L(x=y)  > 

By  calculating  a  (=  a)  and  using  Fig.  59,  [12]  we  at  once  obtain  good 


(3.18) 


estimates  of  s. 


(=  Q).  For  the  large  rms  underice  profile  slopes  en- 


c  j  u  i  mu  voj  i  -  bdC  k  '  ^  '  -  °  -  - -  ■  —  —  r  •  ~  ■  ■  ■  —  -  ■  ~  r 

countered  here  we  see  that<^ac)>can  be  noticeably  reduced  from  unity 

(no  effective  shadowing),  vide  Table  3.2  (and  Fig.  3.1),  again  unlike 

2 

the  ocean  wave  surface  environments  commonly  encountered,  where  o.  , 


is  comparatively  small,  C,  Sec.  7.4  of  [4],  where  ^ack  ~  1. 


Table  3.2  Backscatter,  9Q j  =  30°;  ( ct>_ -  =  tt/21 


(x*y) 


i  1 

0 

0.176 

0 

3.11  •  10*2 

(-30.4)  9.09  -  10 
(-27.4)  1.82  -  10 

0.364 

0.577 

0.839 

1.19 

1.32  •  10*1 

3.33  •  10'1 
7.04  •  10'1 [ 

1.42 

(-23.2)  4.77  -  10 
(-19.7)  1.06-10 

(-16.7)  2.15-10 
:  (-13.7)  4.24-10 

4.20  •  10  - 

...LUIL. 
0.826 
1.37 
1  95 
2.35 


(-30.4)  9.09  -  10" 


(-13.7)  4.27*  lO'^j  (4.84*  10,  0.80 

(-8.2)  1.50  •  10'1  |  13.06  •  10* *|  0.63 

(-4.3)  3.72  •  10'1  j  12.10  •  10* M  0.49 

(-0.7)  8.52  •  10*1  '  il .48  •  10*H  0.37 
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■ 


.* 


For  the  backscatter  condition  taken  here  (0  j  =  30°;  <t>0j  =  11 / 2 ), 
Eqs.  (3.12b),  (3.16)  reduce,  cf.  (A. 2-13),  to 


N 


(0) 

LH-inc 


0ot=8oo; 
(•••  <*>=10° ) 


N 


16 


(0) 

LH-coh 


16 


back 

♦oT=”/2 


back 

vr/2 


cos4  0qT  (cr2tan20oT  +  l) 

9.09  •  10"4  {(32.16)a2  +  l)J  , 


-8cos6oy{6I^2^  (a2)sin20oT+cos26ol-} 


(2) 


-2i 


1.39  (5.82IV  '+3.02  •  10  *■} 


(3.19a) 


(3.19b) 


Figure  3.1  shows  how  a,  ,  a2,  N.^,.  ,  etc.  vary  with  the  rms  slope  angle. 

L  L  Lii*  i  rlC 

Typical  values  of  <<p£>%  are  0(30°)  for  ice  bottom  profiles,  [16].  As 

•y  1 

expected,  increasing  <4>  >2  leads  to  increasing  tilt-factors  and  smaller 

o 

shadowing  coefficients  (S  ),  because  on  the  average,  the  effective  angle 

of  incidence  is  thus  increased  over  9qT  as  a  result  of  the  larger  (bottom) 

slopes.  Also,  for  the  same  reason,  the  fraction  of  the  scattering  sur- 

2 

face  cast  in_shadow  is  increased,  i.e.,  S  is  decreased.  However,  com¬ 
parison  of  S2  (db)  and  N.^.  (db)  from  Table  3.2  shows  that  the  tilt- 

Ln-lnC  «  ^ 

factor  increases  (with  or  <j>  2)  considerably  more  rapidly  than  the 

shadowing  factor  decreases,  so  that  the  net  effect  is  a  significant  (rela¬ 
tive)  gain  in  scattering  strength  due  to  the  increasing  rms  slopes  of  the 
bottom  ice  surface,  as  this  may  occur.  Note,  for  the  ocean  wave  example 
discussed  in  Sec.  3.1  of  [4],  that  the  rms  slope  angle  is  13.1°  (corres¬ 
ponding  to  a  mean  wind  speed  v  =  20  knots),  cf.  Eq.  3.2e-[4],  for  which 

2 

oL  =  0.233,  from  (3.14),  with  5^ac|<  =  -0.4  db,  as  shown  on  Fig.  3.1  here. 
In  this  case  the  associated  tilt-factors  are  comparatively  small. 


3.3  Model  la:  A  Numerical  Example — Incoherent  Backscatter  at 
0oT  =  80°,  fQ  =  1  kHz 

With  the  results  of  Sec.  3.2  in  hand,  let  us  apply  them  to  Eq.  (3.8) 
to  obtain  the  (theoretical)  backscatter  cross  section  for  a  critical 
geometry,  which  tests  the  adequacy  of  our  physical  model.  Such  a  test 
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■3*1  J" 


TD  7375 


occurs  at  "high  frequencies"  and  small  grazing  angles,  where  the  geo¬ 
metrical  acoustics  (i.e.,  "facet")  component  of  (3.8)  is  negligible  vis- 
a-vis  the  Bragg  term,  [4].  Accordingly,  we  choose  0oy  =  80°  (<j>  =  10°) 
with  4>oT  =  n/2,  cf.  Fig.  A. 2-la,  and  signal  frequency  fQ  =  1  kHz.  Thus 
we  have  the  following  calculations: 

A.  Calculations  (0  y  =  80°;  f  =  1  kHz;  <4>2>^  =  0°  -  50°) 


(1.)  kQ  =  2it/X0  =  2irf0/c0  =  2tt  103/1.5  •  103  =  4.188  (rad  nf1) 
•••  k4  =  3.08  •  102  (rad  rrf1)4 


(3.20a) 


(2.)  9nT  =  80°:  coseoT  =  0.174;  cos40oT  =  9.09  •  10-4; 

taneoT  =  5.67  ;  tan20QT  =  3.21  •  10  ; 

•••  2k  sin0„T  =  8.25 
o  o  I 

(3.)  From  (3.10)  and  Table  3.2  we  qet  fi?  S2  for  <4?>'2  = 

O  L. 


sin0Qj  =  0.985  \ 
sin40oT  =  0.941 j 

(3.20b) 

0,  10°,  ...,  50°. 

(3.20c) 


(4.)  WH(0,2kQ  sin0oT),  (2.11b): 
4-rrC 


WH  =  r,2 


Ckc+(2kosin0oT)"] 


1  —  2  2  ’  kc  =  0,05  m  ’ 


C1  =  0.013 


TTC, 


4k4sin40  , 
o  oT 


(1.09  •  10 


)• 


(3.20d) 


B.  The  Backscatter  Cross  Section, 9  ^cl; 

Combining  the  numerical  results  of _A  above  with  those  in  Table  3.2 
and  applying  these,  in  turn,  to  Eq.  (3.8)  for  the  backscatter  cross  sec¬ 
tion,  specifically  at  small  grazing  angles  (0qT  =  80°),  gives  us  the  re¬ 
sults  shown  in  Table  3.3  following: 
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Table  3.3  Model  Ia--Backscatter  Cross  Section, 
- 1  nc_ 

[  eoT  =  80°;  (4>qT  =  V2);  fQ  =  1  kHz] 


Eq.  (3.14)  Eq.  (3.8) 


Eq.  (3.8)  (no  Kirchoff  correction) 


As  expected,  the  frequency-independent  "facet"  component,  of 

(3.8),  is  quite  negligible  vis-a-vis  the  Bragg  component,  u, 
until  the  rms  slopes  become  large  <<P^>^  ^  45°(e.g.,  ^  4  1),  at  1  kHz 
specifically.  For  example,  at  =  1.42  (C^*5  =  50°)  the  facet  com¬ 
ponent  actually  dominates  as  shown  in  Fig.  4.1.  Also,  because  of  the 
general  form  of  the  wave  number  spectrum  (2.11b),  the  Bragg  component 

is  essentially  independent  of  frequency,  cf.  (3.20d)  in  (3.8) 

•  ri  c  *  n  90 

(provided  the  effective  signal  frequency  fQ  is  such  that  (2kosin0oy)  >>k£ 

=  (0.05)^,  cf.  (3 . 20d ) ) .  Then,  4  4iTk~4C^/sin^0oT.  Here,  at  the  smaller 

grazing  angles  (0nj>7Oo)  this  means  f  0(>  100  Hz).  The  backscatter 
(0)  0 

cross  section  a : n '  for  Model  la  is  shown  in  Figure  4.1  ff.  The  curves 
x  x-  x  x  include  a  +2.5  dB  Kirchoff  correction,  taken  for  0qT  =  80° 
from  Fig.  13  of  [19].  This  is  appropriate  for  rms  (wave)  slopes  of 
13.1°  [=  =  5.42  10"^,  for  mean  wind  speeds  of  v  =  20  knots  (=  10  m/sec)], 

but  is  probably  too  large  for  the  larger  ice  slopes  here,  since  the  larger 
"tilt-factors"  reduce  the  penumbra  effects  of  near-grazing  incidence. 

Thus,  the  corrected  curve  for  nc -back  1S  a  conservatlve  tipper  limit. 


C.  Remarks  and  Conclusions 

Our  results  are  not  sensitive  to  the  "splitting"  frequency  f^  (or 
wave  number,  k^  =  2^/^q),  cf.  [13].  However,  it  is  reasonable  to  choose 
this  frequency  less  than  the  signal  frequency  ( f Q ) ,  in  any  case  such  that 
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Figure  4.1  Backscatter  Cross  Sections  at  9ox  =  80°;  f0  =  1  kHz:  Model  lb: 

Eq.  (4.5),  Model  la:  *-x-x-,  Eq.  (3.8):6(0)_h  ±  3-^. 

(All  have  Kirchoff  correction  of  +2.5  dB,  except  - ,  Eq.  (3.8).) 
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2  2  2 

the  Rayleigh  number,  RH  (=  aHbQk0)>  for  the  reSL,'ltin9  "high-frequency" 
Bragg  component  is  small  compared  to  unity,  so  that  the  associated  per- 
turbational  expansion  of  the  Kirchoff  model  (in  the  manner  of  [4])  is 
valid  for  obtaining  this  Bragg  contribution.  [See  the  remarks  at  the 
beginning  of  Sec.  3.1.1 

The  wave  number  spectrum  used  for  splitting  here  is  given  by  (2.11b). 

p 

We  calculate  cr^j,  the  mean  intensity  of  the  Bragg  component,  from 


2tt 


°H  =  /  WH(k)dk/(2Tr)^  =  /  d<j> 


4ttCx 

(kV)2 


kdk 

(2tt)2 


=  a2  [l  +  (kD/kc)2]"1 ;  a2  =  (yk2  =  5.68  m2 ;  kc  =  0.05  rad  m"1  . 

(3.21) 

2 

Table  3.4  shows  Ru  =  (aub  k  )  for  various  “splitting"  wave  numbers,  kn, 

H  H  O  O  U 

at  a  signal  frequency  of  1  kHz  and  10°  grazing  angle,  in  the  backscatter 
regime,  for  the  Lyon  data  of  Sec.  2. 


Table  3.4  Model  la:  Rn  g  f&  =  1  kHz,  9qJ  =  80°  (y  =  tt/2) 


fD  (Hz) 

XD  (m) 

H 

!  30 

0.126 

Rh  =  1.63 

50 

30 

0.209 

■ 

0.650 

100 

15 

0.419 

fl  B 

0.619 

;  200 

7.5 

0.838 

2.02  •  10"2 

4.27  .  10"2 

400 

i  3.75 

1.68 

5.03  . 10"3 

1.06  .  10"2 

|  1  kHz 

1.5 

4.19 

_ _ _ _ _ 

1 

o 

• 

r- 

o 

00 

1.71  .  10"3 

2  kHz 

0.75 

-  ,  —  ■ 

8.38 

— 

2.02  •  10'4 

4.27  •  10"4 

f 

For  the  above  geometry  (and  data)  choosing  fg  >  0(100  Hz)  yields  a  satis¬ 
factorily  small  Rayleigh  number  (fn  ^f  -)  for  the  Bragg  component.  In 

u  nl 

the  present  situation,  i.e.  (3.8)  both  tbe  "^3cet"  and  Bragg 

components  are  independent  of  kg  (>  kQ). 


29 


TD  7375 


Conclusion:  It  is  quite  apparent  from  Fig  4.1  and  Table  3.1  that 

the  simple,  single-sur^ jce  model  (Model  la),  however  arti¬ 
ficially  split  to  gi  'e  "facet"  and  Bragg-scatter  contribu¬ 
tions,  yields  too  small  (incoherent)  backscatter  cross 
sections,  a  these  comparatively  small  grazing  angles 

Land  intermediate  and  "high"  frequencies,  f  >0(100  Hz)]. 

Even  allowing  for  the  large  experimental  spreads  (a  =  0(15  db)) 
in  the  data  (Fig.  4  of  [  5]  ,  and  Milne  [20],  Mellen  and  Marsh 

[2lJ),  shown  as  the  dotted  line  in  Fig.  4.1  here  for  =  80°, 

f  =  1  kHz,  Model  la  is  inadequate  to  account  for  the  ob¬ 
servations. 

Accordingly,  we  must  investigate  the  more  complex, 
two-component  surface  Model  lb,  described  in  Sec.  2.2  above, 
to  see  whether  it  can  reasonably  account  for  the  larger  ob¬ 
scured  backscatter  intensities. 

Finally,  we  note  that  major  numerical  differences  with  the  analogous  ocean 
wave  surface  models  stem  not  from  the  analytical  results,  cf.  (3.8),  but 
from  the  much  larger  rms  slopes  encountered  in  the  ice  profiles,  as  well 

as  from  the  empirical  details  of  the  profile  data  themselves  (cf.  Secs. 

1,  2). 

4.  Scattering  Cross  Sections  for  Model  lb 

As  we  have  seen  in  Section  3,  a  single-component,  or  single-surface 
model  appears  to  be  inadequate  to  explain  tne  observed  (back)scatter 
cross  sections,  cf.  Fig.  4.1.  This  strongly  suggests  that  an  additional, 
but  independent,  small-scale  surface  structure,  riding  on  the  original, 
larger-scale  underice  surface  may  provide  the  required  scattering  mech¬ 
anism,  which  accounts  for  the  empirical  scatter  cross  sections,  particu¬ 
larly  at  the  smaller  grazing  angles  (6^  >70°). 

4.1  Scattering  Cross  Sections;  Backscatter  Cases,  Model  lb 
Thus,  we  modify  Model  la,  (3.5),  according  to 

=  Iz(q_  +  ?H*  +^Ach  (  =  Model  Ia  +  il ACh ^  ’  t4-1) 


it>  *  +£ A^h 
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where  the  additional,  smal  1  -scale-surface  component  is  ^  =  £h(_r),  which 
is  now  regarded  as  normally  distributed,  =  0,  vide  Sec.  4.3  ff. 

The  unit  normal  _n  ,  cf.  (2.15),  reduces  to 


((i  C  +  ?  4  -  i  )  / /l  +  } 

X  «*y  y  f  X  y 


Cx  =  3C/3x,  (4.2) 


the  contribution  of  the  "high-wave  number"  component,  4^,  being  ignor- 
able  vis-a-vis  that  of  in  4^.  [ T r.  fact,  we  "ecognize  (4.1),  with 
_Qa  *  ,  as  being  formally  analogous  to  our  previous  ocean  wave  surface 

model  (2.1)-[4],  which  includes  a  soliton  or  haudraulic  jump  component 
riding  on  the  underlying  gravi ty-capi 11  ary  wind-wave  surface.  Here, 
however,  may  be  "two-sided,"  i.e.,  =  o5  unlike  the  soliton  model 

used  in  the  wind-wave  surface  study  [4],  cf.  (4.11),  (4.12)  ff.] 

Our  analysis  here  consists  of  a  slightly  modified  version  of  the 
treatment  described  in  [4].  The  incoherent  scattering  cross  section 
may  be  compactly  given  by 


i  ?(0)  s  {3(0)  +  3(0)  }  +  3(0) 

inc  inc-L  inc-H  inc-h  ’ 


(4-3) 


where  the  first  two  terms  of  (4.3)  are  respectively  provided  by  (3.6) 
above,  in  the  general  bistatic  cases,  specifically  under  the  Kirchoff 
conditions  described  at  the  beginning  of  Section  3,  which  for  these 
classes  of  profile  data  are  applicable  for  illuminating  frequencies 

f  =  0(2  50  Hz),  viz.,  the  "intermediate"  and  high  frequencies.  Spe- 

°  .  ,,  *(0)  ... 

cifically,  c71-n^.h  is  given  by 


4 

3(0)  =^2^2_Vn(0) 

inc-h  0  ...  2  lh 

16tt 


(201 


).  W.  (2a 

'  i  nr  h ' _ _ r 


■—0  |x,y'inc  hv-~o|x,y 


,) 


(4.4) 


cf.  (2.13)-[4],  where  is  the  wave  number  intensity  spectrum  of  the 
additional,  independent  surface  component, 

For  backscatter  we  use  (3.8)  for  the  first  two  terms  of  (4.3)  and 
write  (4.4),  4>oT  =  */2,  as 
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To  obtain  an  estimate  of  the  wave  number,  k<.,  at  which  a  new  spectrum  is 
required,  we  solve  (A. 1-13)  or  (A .  1-  13a )  for  k<-,  getting  (A. 1-14),  >>1. 

The  results  are  tabulated  in  Table  A. 1-1: 


Table  A. 1-1  Spectrum  Transition  Wave  Number 
k$,  C1  =  0.013;  kc  =  0.05  (m_1) 


I  2  lx 

1  <4>|>  2:  rms  slope 

H1BS39 

k$  (rrf1) 

As(=  2ir/ks)  (m) 

0° 

0 

0 

10 

3.11  * 10'2 

±  0.47 

±  13.3 

20 

-1 

1.32  *  10  * 

2.18  -103 

2.88  • 10~3 

30 

3.33  -10"1 

1.10  •  1010 

-10 

5.71  • 10  iU 

40 

7.04  * 10_1 

50 

1.42 

_ 

_ 

Clearly,  the  transition  wave  number  k<-  is  much  too  large  for  the  rms 
slopes  encountered  in  practice  with  ice  profiles,  e.g.,  0(30°)  so  that 
k^  =  O(1010)  or  =  O(l-rl0'10  m).  This  supports  the  argument,  developed 
in  the  text  (esp.  Secs.  3,  4)  that  a  true  two-component  or  two-scaled  pro¬ 
file  surface  should  be  considered,  since  a  single-scaled  profile,  however 
arti f icial  ly  split,  yields  much  too  small  backscatter  cross  sections  (for, 
in  the  text,  grazing  regimes).  As  we  note  above,  such  a  profile  requires 
much  too  short  a  cut-off  wavelength,  a  to  provide  the  magnitude  of  rms 
slopes  typically  encountered,  vide  Table  A. 1-1  above. 


A. 1-4  Covariance  Functions  and  Wave-Number  Spectra  (Model  lb) 

For  our  surface  elevation  model  (4.10),  when  each  projection  element 
is  independent,  so  that  a  poisson  model  of  element  location  may  be  ap¬ 
plied,  [20]- [22] ,  we  readily  find  that  the  covariance  of  the  profile  is 
given  by 


Khdr) 


/  P(r ' 
A  ~ 


)  <c(irL.'  ’®)^([i+lr"ri’i)>  ir' 

0 


(A. 1-15) 
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since 


jitl,  fj  Hk JS  fk)  , 


where  is  the  amplitude  wave-number  spectrum  of  C(r).  Thus,  we  write 


•"  (If) 


and  since 


\(x=y)  /  $  kxlkx2 


-iJ<1*r1+'fJ^2*XL2  & 


i^1 
*)2 1 


(A. MO) 


S^)S^F=  (2TT)^(k2-k1)W2(k1)?,  with  Ar  =  rz-rv  (A. 1-11) 

2  2  2 

which  defines  the  wave  number  spectrum  intensity  W^,  we  use  k  =  k„+k.. 

p  c  a  y 

=  2V.  and  (A. 1-11)  directly  to  obtain  the  result  (A. 1-9).  Equivalently, 
x  y 

for  these  isotropic  cases 


K(x=y)  ’  (f )  *  (|)  '  IJ  ^2  =  h  /’  kH<*>- 


(A. 1-12) 


Now,  using  (A. 1-7)  in  (A. 1-9)  and  truncating  the  spectrum  Wj  at  some 
(large)  spectrum  transition  wavenumber  k<.,  to  get  a  finite  rms  slope,  we 


have  here 


L(x=y) 


=  Cj  . =  T  !log(KSc+1)  +  ^2~  ’  1 

0  (k2+k2)2  2  (  5c  1+kL 


KSc=  kc» 


(A. 1-13) 


T{2  log  KSc  -  15  ’  he”  l' 


(A. 1 - 13a ) 


kS  ’  kc  expl5-  [1  *  2aL(x-y)/Cl])  • 


(A. 1-14) 
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from  Watson  [10],  p.  410,  (2),  with 

W2(k)Menen  =  2C1/(kj;+k2)2  ,  cf.  (4)  of  [5].  (A.l-7a) 

Similarly,  using  (A. 1-7)  in  (A.l-5a)  gives  the  correct  covariance, 

viz. : 


°cpcte> 


rr/  / 


4itC 


1 _ i  kAr  cos  4) 


(2tt)2  -J  (k2+k2)2 


dkxdky 


Zn  00  -i  *  ,  4TrGn 

ikAr  cos<t>  _ 1 


-^—5-  /  d<p  f  kdke1 


(2^)  0  0 


(k2+k2)2 


00  J  (kAr)kdk  /C. 

=  2C  /  ■■  2-2T=  TlkcMl(kcAr) 

1  0  (k2+k2)2  \k2  '  c  1  c 


(A. 1-8) 


as  required,  cf.  Watson  [10],  p.  434,  (2). 


A. 1-3  RMS  Slopes  (Model  la) 

The  calculation  of  rms  slopes  o 


L(x=y) 


of  the  ice  profile  can 


readily  be  carried  out  once  the  spectrum  (or  the  correlation  function) 

is  known.  From  (A. 1-4),  or  (A. 1-8),  it  is  at  once  evident  that  with 

2 

this  spectrum  model  aL(x=y)  becomes  infinite,  so  that  for  suitably  small 
wave  numbers  the  empirical  spectrum  clearly  is  inadequate,  i.e.,  does 
not  fall  off  fast  enough  as  k  -’-00.  Thus,  we  may  use  the  empirical  spec¬ 
trum  (A. 1-7)  up  to  some  effective  cut-off  frequency  k^,  at  which  point 
another  spectral  law  takes  over.  Just  what  the  new  spectrum  (at  k>  k^) 
is  remains  to  be  determined  by  experiment  [cf.  our  remarks  a  propos  of 
Model  lb.  Sec.  4]. 

For  the  (isotropic)  slopes  we  have,  in  general, 


(A. 1-9) 


43 


and  since 


1  im 


g  xK^(x)  =  1,  we  have  at  once  (cf.  Watson,  [10],  p.  185): 


4  -  ci/kc 


=  0. 013/(0. 05)2  =  5.20  (m2) ;  *‘-c^=2.28[-2. 38=6^  ,  (2.6)]' 


(A. 1-4) 


;  P;(Ar)  =  I^ArK^k^r), 

< _ _ _ I _ 1 

which  is  Eq.  (2.10b). 

A. 1-2  (Isotropic)  Two-Dimensional  Cases  (Model  la) 

For  the  associated  spectra  and  correlation  function  in  the  two- 
dimensional  cases  we  have,  generally. 


(2-dimensions):  K_(Ar)  =  a2o  (Ar)  =  f  f  W„(k)e"ik’~  , 

“  ?  5  ~  -co  c  ~  (  2tt  )  ^ 


with 


W?(k)  =  I  I  K  (Ar)elk‘-  d(Ar) 


Similar  to  (A. 1-2)  there  is  the  relation 


^V^Mellen.tS]  " 


(A. l-5a ) 


(A. l-5b) 


(A. 1-6) 


Let  us  use  (A. 1-3)  and  the  postulated  i sotropic  character  of  the 
wave  number  spectrum  to  determine  the  two-dimensional  form  of  W^^)  here, 
viz.  (2.11b).  From  (A.l-5b)  we  have 

2  00 

w„(k)  =  /  dd>  /  Ar  d(Ar)a2(k  Ar)K. (kAr)eiAr*  C0S(p 
0  0  C  c  i  c 

=  2no£  f  ArJo(kAr)K1(kcAr)(k(.Ar)d(Ar) 

oo 

=  (2mj2/k^)  /  JQ(kx)K1(x)x2dx,  b  e  k/kc , 


,  0^4>^  2;  k  =  /k2+k2  ,  Eq.  (2.11b),  (A. 1-7) 

a  y 
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Appe  n  d  i  x__Ai  1 ...  _  Jja  v  e  Jlunibgji^e^Xra_^and .  S  qa  t  i  a  1 
Correlation  Functions.  Models, la,  lb 


A. 1-1  One-Dimensional  Correlation  Functions 
and  Spectra  (Model  la) 

For  the  wave  number  spectrum  and  associated  spatial  correlation 
function  we  employ  the  following  Fourier  transform  relations: 

oo 

(1-dimension) :  l^(|Ar|  )  =  a2p^  (j Ar|  )  =  /  W1(k)e'ikAr  ||,  Ar  =  |Ar|  , 

(A. 1-la) 

with 

VUk)  =  J  k  (Ar)elkArd(Ar),  (A. 1-lb) 

1  -  oo  ^ 

where 

CO 

=  K^(0)  =  /  Wx(k)  |jf-  (A.l-lc) 

for  the  mean  intensity  of  the  profile  elevation,  C .  In  reference  [5] 
the  wave  number  intensity  spectrum  is  defined  as 

2"Wl<k>Mellen,[5]  *  “l<k> 


here. 

It  is  instructive  to  obtain  K^(Ar)  from  the  empirical  fit  of  Mellen 
et  al.  [5],  (2.9)  viz. 


Wj(k)  =  TrC1/(k^  +  k2)3/Z  . 


(A. 1-3) 


Thus,  using  (A. 1-3)  in  (A. 1-la)  gives* 
*C 


'1 


afp r(Ar)  «  /  o  ,  V,  t 

■  ‘  (k2+k2)3^2  ^ 


dk  -ikAr  _  r  r 

“  e  -  C.J 


/  C°S  ^3/?dk  =  4(k  Ar)K  (k  Ar), 

1  0  (k2+k")3/2  k2  c  1  c 


(A. l-3a) 


*In  equation  (2)  of  [5]  "constant"  should  be  replaced  by  "(1/2) 
constant."  This  also  insures  that  (4)  of  [5]  agrees  with  (2.11b),  etc. 
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(i).  To  include  specific  keel  structures  in  the  manner  of  Models 
IIa,b  (Sec.  2.2); 

(ii).  To  compare  the  results  of  (i)  with  those  of  the  present 
study; 

(iii).  To  evaluate  forward,  as  well  as  backscatter,  including  co¬ 
herent  scatter  components. 

Critically  needed  now  are  active  experimental  data,  both  linear  and  di¬ 
rectional,  at  wavelengths  0(6  m  -  1  cm)  on  underice  profiles  for  typical 
regions  in  the  ocean.  With  such  data  we  can  evaluate  and  compare  the 
various  scattering  models,  and  select  the  one  appropriate  to  the  profile 
regime  in  question. 
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present  results  depend  primarily  on  the  chosen  parameter  values  of  the 
model  in  question, and  only  secondarily  on  the  specific  analytic  structure. 

The  conditions  under  which  our  quantitative  models  above  have  been 
constructed  and  evaluated  are,  concisely: 

Model  Conditions 

(1) .  Far-field  operations,  for  both  source  and  receiver  vis-a-vis 

the  underice  bottom  (Appendix  A. 2); 

(2) .  Kirchoff-Bragg  scatter  conditions  (Sec.  3),  requiring  "inter¬ 

mediate11  0(0.10-1.0  kHz),  and  "high"  frequencies  0(5*  1.0  kHz), 
which  are  in  the  large  Rayleigh  number  regimes; 

(3) .  Composite  gaussian  surface  profile  model  (Secs.  2,  2.1),  based 

on  Lyon's  data  [ l]  (but  see  Sec.  2.1  for  limitations); 

(4) .  An  isotropic  covariance  (and  associated  wave-number  spectra) 

for  the  large-scale  profile  component,  based  on  Mellen  et 
al . ‘s  recent  work  [ 5]  ; 

(5) .  Finite  (rms)  slopes— Sec.  3 . 3  C_;  gaussian  pdf  of  slopes; 

(6) .  "Single-bounce"  geometry  for  ice,  source,  and  receiver;  and 
'  Vc  «  0; 

(7) .  "Hard"  ice  boundary,  with  local  plane  wave  reflection  coef¬ 

ficient  (i?0=l)»  cf.  Sec.  3.2,  A. 

(8) .  Large  "tilt-factors"  and  non-negl igible  shadowing  corrections. 

Sec.  3.2,  4.2. 

(For  details,  see  the  text  above  and  Appendixes  A.l,  A.2).  Unlike  the 
related  problems  of  scattering  from  the  ocean  wave  surface  [4],  [12], 
the  surface  slopes  involved  here  are  much  larger,  0(2. 5-3),  with  conse¬ 
quent  modification  of  the  analytical  and  numerical  results,  cf.  Secs. 

3.2  B,  C;  4.2,  Table  3.2  above. 

Our  Models  la,  lb  do  not  explicitly  designate  a  keel  structure, 
for  the  reasons  we  have  already  noted  in  Section  1.  Neither  are  we  con¬ 
cerned  here  with  "low"-frequency  scatter  [f  =  0(«S  30-50  Hz)],  and 
multiple  surface  interactions,  produced  by  suitable  geometries  in  a  non¬ 
zero  gradient  ocean  (Vc  f  0).  These,  and  other  model  situations, are 
reserved  to  subsequent  studies.  Of  these,  our  immediate  next  steps  are 
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is  not  critical;  rather,  it  is  the  values  of  height  (~o  )  and  scale 
(~&h)  which  are  significant  (except,  of  course,  at  the  very  high  wave 
numbers,  which  influence  the  (rms)  slopes). 

5.  Concluding  Remarks 

The  principal  general  results  of  this  study  are: 

I.  Determination  of  the  incoherent  (backscattering)  cross  sec¬ 
tion,  for  Model  la:  a  single  (isotropic)  underice  sur¬ 

face  (Sec.  2,  Sec.  3); 

II.  Extension  of  Model  la  to  Model  lb,  which  adds  an  additional 
surface  component,  of  much  smaller  scale. 

The  general  result  II  is  generated  from  the  observation  that 

III.  Model  la  is  shown  to  be  inadequate  to  account  for  the  empi¬ 
rically  observed  backscatter  cross  sections  (Fig.  4.1), 
particularly  in  the  critical  regimes  of  small  grazing  angles 
(0ot>7o°)  and  "high"  frequencies  fQ>0(o.5  kHz).  Further¬ 
more,  Model  la,  if  applied  for  all  wave  numbers,  yields  in¬ 
finite  (rms)  slopes  (Sec.  3.3 _C). 

IV.  Model  lb  can  account  for  the  observed  (back)scatter  data, 
cf.  Fig.  4.1.  The  scale  of  the  added  independent  surface 
component,  appears  to  be:  rms  elevation  ah  =  0(8-  15  cm); 
e'^-correlation  distance  =  0(4-7  cm).  (Model  lb  always 
has  finite  rms  slopes.) 

V.  These  models  and  results  strongly  support  the  need  for  pro¬ 
file  data  at  wavelengths  0(6m  -  1  cm),  from  which  appropriate 
wave  number  spectra  can  be  constructed  for  model  use. 
Accordingly,  a  two-,  independent-component  model  of  the  (underice)  pro¬ 
file,  involving  small-  and  large-scale  surface  elements,  can  account 
for  the  observed  (back)scatter  cross  sections,  particularly  in  the 
limiting  regimes  of  "high"  frequencies  and  small  grazing  angles.  Whether 
such  a  model  (Model  lb  here)  is  appropriate,  i.e.,  adequately  repre¬ 
sents  the  underice  scattering  mechanisms  involved,  remains  to  be  de¬ 
termined,  both  theoretically  and  empirically  [the  latter  from  V  above, 
and  remarks  below:  "next  steps"].  One  limitation  here  is  the  restriction 
of  ;A,  (4.1),  to  a  single  gauss  process,  albeit  with  the  correct  (i.e., 
empirical)  second  moments,  as  described  in  Sec.  2.1  above.  Models  II 
[Secs.  2. 2-2. 4]  are  not  so  constrained.  It  is  emphasized,  also,  that  the 
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-x..  _ 

/  *V*c 

\(  e”*a. 

n 


Ar  -> 


Figure  4.2  Sketch  of  K.  (Ar),  Eqs.  (4.12),  (4.12a);  l  =  "correlation 

distance,"  K,  =  e~  a.  . 

h  h 


Remarks 

Figure  4.1  shows  the  incoherent  backscatter  cross  section  .  , 

mc-h 

attributable  to  the  small-scale  component  riding  on  the  large-scale 
profile,  cs,  for  a  selection  of  parameter  values  (oh,£h).  We  have 
chosen  values  which  give  cross  sections  in  the  range  of  empirical  values 
when  the  rms  slope  is  0(30°).  Typical  values  arech  *  0(10  cm)  and 
£h  =  0(4-5  cm),  or  £ c  =  (5+  -  7  cm),  for  our  present  model.  Here 

^inc-back  ~  ainc-h  for  s^opes  °(20c>-  45° ) :  the  independent,  "small- 
scale"  component  riding  on  the  large-scale  profile  is  here  clearly  the 
dominant  contributor  to  the  incoherent  backscatter  cross  section,  for 
these  (comparatively)  small  grazing  angles  and  high  frequencies, 

D ( f Q  =  1  kHz).  The  rms  roughness  of  ch,  as  measured  by  oh,  is  0(2-5%) 
of  that  of  the  (composite)  large-scale  profile,  [Of  course,  at  this 
stage  we  do  not  have  experimental  verifications  of  the  scale  of  these 
numbers,  except  indirectly  through  the  postulate  of  Model  lb  and  the 
comparisons  shown  in  Fig.  4.1.] 

Finally,  if  we  use  the  generalized  version,  (A. 1-22),  of  the  wave- 
number  spectrum  noted  above  (4.14),  we  obtain  for  the  same  parameters 
(and  d z  =  /it/ 2  £^)  essentially  the  same  values  of  ^ j nc-h *  more  Precisely» 
values  only  0(2. 5-3.0  db)  lower  than  those  shown  in  Fig.  4.1.  This  is 
in  keeping  with  our  remark  above  (beginning  of  Sec.  4.3)  that  the  "shape" 
of  the  individual  projections  forming  the  small-scale  surface  component 
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4.4  Model  lb.  A  Numerical  Example: 

Backscatter  g  3_T=80  t  <|).t=tt/2;  f_=l  kHz 


Backscatter  @  3_T=80  , 

- OT - =— 7T 

Combining  (4.7)  for  N,\ 


_  Combining  {4_.7)  for  N^;  cf.  Table  3.2,  (4.14)  for  W^(O,2kosin0oj) ; 

,  (3.10)  and  S  ,  Table  3.2,  in  (4.5)  gives  us  the  desired  quantitative 
0  (0) 

structure  for  the  backscatter  cross  section  d)  1  .  ,  arising  from  the 

mc-h  3  /0\ 

added  surface  mechanisms,  c  .  Table  4.1  shows  a  range  of  values  of  8)  '  ,  . 

h  3  mc-h 

for  various  choices  of  model  parameters 


Table  4.1  Model  Ib--Backscatter  Cross  Section 
^inc-h’  (4.5) :  (Uncorrected*) 
[eoT  =  80°;  fQ  =  1  kHz;  (4>oT  =  ir/2)] 


. 

Uh  •  3.0  cm 

r 

rms  slope 

|oh  •  3.0  cm 

6. 

|  -68.6  (db) 
|  -58.4 


10.0  !  15.0 


-5 


•we  add  a  nominal  +2.5  db  to  each  entry,  to  provide  the  needed  correction  for  the 
at  coT  •  80°  (e*10°),  from  Fig.  13  of  [19]  ;  cf.  Fig.  4.1. 


approximation 


ul 


The  shaded  regions  of  the  Table  indicate  cross-sectional  values  that 
fall  within  the  A-region  (±15  db)  of  uncertainty  in  the  empirical  data 
shown  in  Fig.  4.1.  Also  in  Fig.  4.1  are  four  curves,  representative 
of  possible  model  values  (Model  lb).  (The  A-region  is  -17  to  -47  db.) 

Figure  4.2  shows  K^(Ar)  and  typical  parameter  relations:  in  our 
present  model  lQ,  and  therefore  l  are  fixed;  (for  a  model  in  which 
Lq  (*  22.^)  is  random,  see  Appendix  A. 1-4).  Here  i  (=  /2£h)  is  the 
correlation  distance  of  the  profile  component  (;.  ). 
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cf.  Appendix  A. 1-4,  with 

a2  =  Kh(0)  =  Y0  L2tt/8;  lh  =  lo/2  ;  Ar2  =  (Ax)2  +  (Ay)2,  (4.12a) 
where  the  mean  value  is  explicitly 

«h>^o'^/4.  (4.12b) 


(Note  that  if  5  =  0:  the  mean  "bump"  height  is  zero,  then  there  is  no 

"dc"-component  in  the  wave  number  spectrum,  cf.  (4.15)  ff.)  Equation 
(4.12)  is  shown  in  Fig.  4.2.  (However,  for  an  alternative  model,  extend¬ 
ing  (4.12),  see  Cases  I, II  in  Appendix  A. 1-4,  and  closing  remarks  in  this 
Section. ) 

The  two-dimensional  wave  number  spectrum  corresponding  to  (4.12)  is 
obtained  from 


Wh(k)  =  /  /  Kh(Ar)e1-*~  d(Ar) . 


(4.13) 


cf.  (A.l-5b).  We  have,  for  £  *  (0,2kosin6oy)  specifically  here 

V*>  =-2™^e'2<k°V1"9°T>2 


(4.14) 


for  the  "continuum"  part  of  the  spectrum.  If  CQ>  0,  then  the  complete 
wave  number  spectrum  is 


=  2  7-2  tt 


“h<i>  *  *0  C  V  Lo^-°>  +  2"°h*he 


•2(k0ihsineoT)' 


(4.15) 


08  ik  Ax 

where  6(k-0)  =  6 (k  -0)6 ( k  -0)  and  /  Ae  d(Ax)  =  2rrA<5(k  -0),  etc. 
“  x  y  x 


(Note,  incidentally,  that  1 ;h  has  finite  rms  slopes,  unlike  above; 
(Appendix  A.l).  The  "dc"-component  in  (4.15)  does  not,  of  course,  in¬ 
fluence  the  scattered  radiation. 
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given  by  the  results  of  Table  3.2  and  Fig.  3.1,  cf.  Sec.  3.2,  C..  The 
dominant  effect,  however,  here  is  that  produced  by  the  "tilt,"  embodied 

NLh-inc>  <4'9>- 

4.3  Model  of  the  Small-Scale,  Additive  Surface  Component,  C. 


A  variety  of  plausible  models  for  the  small-scale  component  is 
certainly  available  at  this  stage  of  our  study.  One  such  model,  which 
we  employ  here,  is  that  involving  random  "bumps"  or  projections,  dis¬ 
tributed  in  space  independently,  riding  on  the  large-scale  surface, 
cf.  (4.1).  These  projections  are  accordingly  poisson  distributed,  and 
may  overlap,  or  ride  upon  each  other  occasionally.  Physically,  these 
"bumps"  represent  inhomogeneities  in  the  local  crystal  formations  on 
the  large-scale  surface,  produced  by  near-surface  freezing  forces. 

A  reasonable  model  of  a  typical  (j— )  projection  is  given  by  the 
gaussian  form 


y.  =  5h(r.rj)  ■  v 


■4|r-ri|2/L^ 


where is  a  (random)  location  (referred  to  the  reference  plane,  SQ, 
cf.  Figs.  A. 2-1,  A. 2-2).  [What  is  primarily  important  in  such  models 
is  the  height  (~C  )  and  spread  (~LQ);  the  actual  spatial  shape  only 
becomes  critical  out  on  the  "tails"  (lx~Xil>>  L  ) •  or  equivalently,  for 
the  high  wave-number  fall  off  of  the  associated  spectrum. J  Accordingly, 
for  the  isotropic  poisson  process  (4.10),  we  may  write  (cf.  Appendix 
A. 1-4)  for  the  mean  and  covariance,  respectively,  generally 

<Ch>  =  V'S^l^  *  M—)  =  Yo<^h^~l^h^l+^‘^>  *  ^4- 

where  tq  is  the  process  density  (in  space)  =  average  number  of  "bumps" 
per  unit  area  (on  CA)  cf.  [20],  [21]. 

For  our  specific  isotropic  model  (4.10),  (4.11)  reduces  to 


Y  ttL2!;2 

S(ir>  ■  V4 


-2Ar2/L2  ,  -Ar2/2f ^ 

6  =  ahe 
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NLh-inc  16  (I<4)sin49oT  +  6l(2)s1n2eoTcos29oT  * 

back  90y=7V2 


(4.7) 


where 


(  «  -(x2+y2)/2a?  }  2 

\  I(0,(a)  a  /  /  S - -  -aVW);  a2  a  vtf?  ; 

\  -*  1  +  +  /  2iraf 

U"y;  (4.8a) 

j  oo  o  2  -(x2+y2)/20x 

J  I(2)(a)  =  /  /  tx  ■  or  y  > - ^ - ^4=  id  ~  I(0)(a)},  cf.  (3.17), 

\  -«  1  +  vr  +  y  2™~ 


(4.8b) 


»  4  4  -(»2+y2)/2o2 

I  ^)(a)  =  /  , - I 

-oo  1  +  x  +  y 


2’°x  a 

■  II7- 1  *  1<0)<a>!- 

a 


(4.8c) 


and  Ei(x)  is  the  (tabulated)  exponential  integral,  [17],  [18].  (Note 
that  with  vanishing  slopes,  a-*-°°,  I^-*l,  and  I^2\  I^+O.)  The  re¬ 
sult  (4.7)  is,  in  the  general  bistatic  cases  (with  4>ny  =  tt/2): 


TT»lfn  *i  “oyl(4)<a>  *  6ao,aozI<2)(a)  +  a42l(0)<a)i-  <4-9) 


cf.  (7.33e)-[4],  and  (A. 2-11),  (A. 2-13,  14)  for  a  a  a  etc. 

/g\  ox  oy  oz 

Values  of  nc/16  for  the  backscatter  cases  are  tabulated  in  Table 
3.2  and  shown  in  Fig.  3.1.  This  tilt-factor  is  seen  to  be  much  larger 
than  NLH-inc^®’  for  the  artific''any  "factored"  surface  =  cL  +  ;H, 
cf.  Sec.  3.1  et  seq.,  as  the  profile  slopes  increase. 

The  other  factors  which  scale  remain  unchanged  in  basic  struc¬ 
ture:  7?  =  0  db,  cf.  (3.10),  and  the  shadowing  factor,  ?,  is  still 
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d(°) 

i nc -h |  back 


■  "l^VtUc  ,  V°-2VineoT> 


R  >>1 
<t>oT=7f/2 


(:,K*  ) 
x>ack' 


As  before,  cf.  (2.14a)-[4],  N^'  is  the  "tilt-factor"  for  associated 
with  the  larger-scale  surface  component  cL»  and  is  given  generally  by 


/n\  /'a  ^  +  a  ?  -  a  )  \ 

n  °>  .  16  /  ox  x  °yy  °z  > 

Lh  inc  \  ^  +  +  2  / 

x  y  L 


so  that  for  backscatter  (with  4>0y  =  tt/2  ) » 

N.^  inr  =  16  i  sin49  ,  ^ - o — j)  +  6sin2Q  TCOS26  T  / - ^ 

Lh-  inc  )  oT\  2  2/  oT  °T  \  2  2  , 

back  l  *x  y  L  x  y  L 


+  cos49 


L ! 


(for  symmetrically  distributed  slopes,  about  <C  >.  =  0),  where  (3.15) 

x  »y  l 

governs  the  distribution  of  C,  £  „  here.  As  for  the  Bragg  component 
(0) 

u  in  (4.3)],  the  pertinent  surface  statistics  for  the  new  com- 
lnc_H  (01 

ponent,  (4.5),  are  the  wave-number  intensity  spectrum:  there  is 

no  requirement  that  in  (4.1)  be  normal,  for  example. 

n  (01 

Finally,  the  complete  expression  for  is  provided  by  (3.6) 
and  (4.4),  generally,  and  (3.8)  and  (4.5)  for  the  backscatter  cases 
treated  quantitatively  here. 


4.2  The  Tilt-Factor  N^h-inc-back 


Again  [cf.  J_,  Sec'.  3.2],  because  of  the  large  rms  slopes  involved, 
we  must  use  the  exact  relations  ( 7 . 64 ) - [4 D ,  (7.65)  -[4]  to  evaluate 
(4.6b),  with  the  help  of  (3.15),  (3.17).  For  the  assumed  isotropic 
surfaces  here  (a,  =  a,  1  (4.6b)  becomes 
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where  .0=  (4  ,  . ..),  i.e.,  all  statistical  parameters  of  4,  and  is  the 
domain  (i.e.,  surface)  over  which  the  individual  projections  can  occur 
(in  our  model,  individual  "bumps"  or  projections  can  overlap).  On  the 
assumption  of  a  uniform  density  of  "bumps,"  (4.10),  on  A,  we  can  rewrite 
(A. 1-15)  as 


Kh(£r)  *  Y0  /  <4(i,0K(A+ir,0)>0  dX 


—  oo  /-4(Xv+^)/Lo-4(A‘+2A  Ax+Aj+2X  Ay+Axfc+Ayfc)/LjV 

=  Vo  /  /  dAxdXy<^  x  y  o  x  x  y  y  V 


(A. 1- 16a ) 

2.„  2w,2. 


L 


(A. l-16b) 


for  these  isotropic  cases.  Evaluating  (A. l-16b )  is  straightforward.  We 
get  here 


~n  /h  -2ir2/L2 

-  V 0<Ct  e 


L?ttx  -2Ar2/L2 


(A. 1-17) 


where  in  our  present  model  we  now  assume  that  w^(l_0)  =  6(Lq-L0),  so  that 
(4.12)  is  the  result. 

However,  it  is  instructive  to  consider  some  other  pdf's  of  LQ: 


Case  I:  W^(L0)  =  Rayleigh  pdf 

Here  LQ,  or  &h,  is  assumed  to  be  governed 


wid0)  = 


Loe 


-Lo /2al 


L0  >  0;  or  w1(4h) 
L. 


by  a  Rayleigh  pdf,  viz . : 

-»f \n°\ 

,  V 

°l  (A. 1-18) 

=  2*h 


Using 

°°  i  -Bx’^-Yx”*3  o  v/2p 

/  xy-1e  dx*f(J)  K  (2/87),  Re(6,y>0)  (A. 1-19) 

n  P  '  ts 1  v 

u  B 

[23],  p.  342,  No.  3.478-(4) ,  we  see  that  (A. 1-17)  becomes,  with  (A. 1-18), 


TD  7375 


Kh(Ar) 


Vo 


l! 


3  ~ll^W/nl 


£he 


dV°£ 


(A. 1-20) 


and  since  ^  |x2<2(x)  =  1,  we  •'ave  directly 


[o\  =  i  -  **) 


(A. l-20a) 


The  associated  (two-dimensional)  wave-number  spectrum  from  (A. 1-5) 
is  now 


W2(k)  =  ^°h°£  x3K2(x)J0(ko^x)dx*  Ar/CTj,  =  x» 


(A. 1-21) 


(A. 1-22) 


by  [10],  p.  410,  Eq.  (2).  Unlike  the  empirical  model  (A. 1-7),  this 
spectrum  supports  a  finite  rms  slope.  From  (A. 1-12)  we  see  that  speci¬ 
fically 


h(x=y) 


_1_ 

4tt 


aZ  4 

/  k3W_(k)dk  = 

0  c 


k3dk 


0  (k2+k2)3 


k2 
2  KZ 
°h  4tt 


°h 


(A. 1-23) 


this  last  from  [23],  p.  293,  3.241-(4)  directly. 

Incidentally,  it  is  not  possible  to  get  the  covariance  (A. 1-3,4) 
from  our  gaussian  model  (4.10),  as  this  would  require  a  pdf  of  &h  of  the 
form 
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W1^V  =  ~  ^2 


,  *h>0. 


(A. 1-24) 


which  is  singular,  i.e.,  not  bounded  or  normal izeable. 


Case  II:  w^(LQ)  =  "one-half"  Gauss  pdf 

In  this  case  L  ,  or  £ h  is  now  assumed  to  obey  a  "half-gauss"  pdf. 


w,U  )  =/^  e  L°/2  *  or  =/^jr  e 

1 '  o '  y  n  o.  1 '  n'  »  7r  a„ 


#z4 


Using  (A. 1-19)  with  (A. 1-25)  in  (A. 1-17)  gives  now 


■  °l  to'-/°l)3/2  S/2  (Ar/°£)- 


Since  ([10],  p.  80) 


K3/2(x)  =/2x  e  [1  +  x  ]  , 


we  note  that 


x3/2k3/2(x)  =  1# 


Consequently,  we  may  rewrite  (A.l-26b)  as 


Kh(ar)  -  (Y0c02  foflj/f  (4f)3/2K3/2(Ar/o,)j 


"  CThkhC&r)  -  o?  &  +  1)  e  4r/<,,! 


where  the  mean  intensity  is  here 


Kh(0)  =  =  (Y0?q  \  a\)' 


JL 


<Lo  =  2V 

(A. 1-25) 


(A. l-26a) 


(A. l-26b) 


(A. 1-27) 


(A. 1-28) 


(A. l-29a) 


(A. l-29b) 


(A. 1-30) 
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as  before,  cf.  (A.l-20a),  and  the  expression,  {  } ,  in  (A.l-29a)  is  the 
normalized  covariance,  k^(Ar). 

The  associated  (two-dimensional)  wave  number  spectrum  follows  from 
(A. 1-29)  in  (A.l-5a).  Taking  advantage  of  the  surface  isotropy  implied 
by  Ar  -*■  |Ar|  =  Ar  here,  using  (A.l-29a)  we  obtain* 

«2(JsJ  -  <4  /"  d(ir)iro2/2  «Wol>e1k4r  C°S*  <*-l-31a) 

//—  oo 

|/Q  J0(ka£y)K3/2(y)y5/2dy  (A.l-31b) 


(k£  =  l/a£) ,  (A.1-31C) 

(A. 1-3  Id ) 


w2(i)  = 


2  2 


2, 5/2  ’  Vk'’ 


n+(ka£)  ] 

erra2k£3/(k2+k2)5/2 


again  from  (2),  p.  410  of  [10].  This  spectrum,  unlike  that  of  the  empirical 
model  (A. 1-7),  also  supports  finite  rms  slopes,  e.g.,  from  (A. 1-12), 


n(x=y)  H  gauss 


=  X  f 

4tt 


k3W2(k)dk 


(I) 


x3dx 


—072 


(=2/3) 


=  ^/°£  (<  “),  (A. 1-32) 

like  Case  I  above.  Here  W,  is  0(k"3)  as  k +«>,  while  W„  D  ,  .  .  is 
_6  2  2-Rayleigh 

°(k  ),  and  W2-empirical  is  °(k  As  yet,  the  various  pdf’s  for  Z^ 

(Cases  I, II)  are  plausible,  but  not  physically  justified  in  any  profound 
way.  One  task  of  a  subsequent  study  will  be  to  attempt  a  rational  choice 
°f  to  be  tested  ultimately  against  experiment.  Basically,  it  is 

the  behavior  as  k-*®  which  appears  to  be  the  controlling  factor  in  such 
model s . 


*We  can,  of  course,  use  (A.l-29b)  directly,  to  obtain  (A.l-31c,d). 
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Appendix  A.2.  Scattering  Cross  Sections: 
Denrntions,  Nomenclature,  and  Conditions, 


There  is  a  number  of  variations  on  the  concept  and  definitions  of 
scattering  cross  section.  It  is  therefore  important  to  define  the 
term  explicitly,  so  that  the  different  definitions  can  be  calibrated 
with  one  another  [which  we  shall  need  to  do  in  order  to  effect  compari¬ 
sons  with  both  the  analytic  results  of  others  and  the  corresponding 
measurements  (cf.  [4],  Secs.  2.A,b  and  Sec.  3  ff.  for  ocean  wave-surface 
scattering)]. 

We  begin  with: 

Incoherent  Scattering  Cross  Sections, 


Here  we  define 


incoh  (scatterin9  at  receiver  R)  .  RQf 


incident 


[at  surface)  A, 


:[(path  loss)  x  (gTgR)2] 


where  1  i ncoh *  incident  ar8,  respectively,  the  intensities  of  the  scat 
tered  and  incident  fields  at  the  points  indicated.  The  basic  concept 
of  the  scattering  cross  section  (for  surfaces)  is  to  eliminate  the 
effects  of  source  level  and  propagation,  i.e.,  the  effects  of  the 
medium--which  are  handled  separately--when  computing  energy  loss,  and 
to  focus  on  the  effects  of  the  random  scattering  surface  itself. 

To  keep  dimensionless,  a  reference  "illumination"  area,  Apr.-. 

is  employed,  whose  specific  form  is  suggested  by  the  composite  beam 
pattern  projection  on  the  reference  or  equilibrium  surface  <c>  *  0:  S 
Figure  A. 2.1  shows  the  relevant  geometry. 


m 


Figure 

The 

( A . 2 - 2a ) 

(A.2-2b) 

(A. 2-2c) 


.2.1  (a)  Backscatter  Geometry  (far-field):  Monostatic  scat¬ 

tering  (ROT),  vide  Fig.  A. 2. 2;  (b)  Bistatic  Scattering 
(R^T):  "forward"  or  oblique  scattering  geometry,  vide 
Fig.  A. 2. 2. 


various  factors  in  Eq.  (A. 2-1)  are  given  by 

incoh  (scat-  at  receiver)  =  =  int;ensity  of  the  (in- 

coherent)  scattered  field  at  the 
receiver  (R)  [cf.  Sec.  7.2  et  seq.; 
Sec.  8  ff.  of  [4]]; 

incident^  surface)  =  K0^°^in/^4lTRoT^2  =  intensity  of  the  in- 

cident  field  at  0^,  on  the  equilibrium 
scattering  surface  SQ; 

gT,gR  =  the  aperture  "gain"  of  the  trans¬ 
mitting  and  receiving  systems,  cf. 
(6.3),  (6.7)  of  [4]; 
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(A.2-2d)  =  a  reference  area  on  the  equi¬ 

librium  surface  SQ,  cf.  Fig.  A. 2.1, 
projected  by  the  composite  T  and  R 
beam  patterns.  (See  A. 2-3  f f . ) ; 

-2aw2c  T 

(A.2-2e)  |  "path  loss"  =  e  0  0  °;  cqT0  =  RqT  +  RqR;  cf.  Sec.  5.3,  Eq -  (5.18)- 
(absorption)  [4];  u)Q  (=2wfQ  )  is  the  (angular) 

frequency  of  the  emitted  signal;  a  is 
an  absorption  coefficient;  cQ  =  (mean) 
wave  front  speed  of  sound  in  the 
water  medium; 

RqR  =  distance  of  the  receiver  (origin)  from  0<., 

Fig.  A. 2.1; 

RqT  =  distance  of  the  transmitter  (origin)  from  0S , 
Fig.  A. 2.1; 

t(r,t)  =  (vector)  wave  surface  elevation,  cf. 

(3.5).  . 

The  definition  (A. 2-1)  for  the  incoherent  scattering  cross  section, 
and  those  similar  to  it  ([3]-[5],  [8],  [9],  [23]-[25],  [29]  in  Ref.  [4] 
here,  for  example)  are  formally  independent  of  range  (R  _,  R  ).  This  is 
not  an  inherent  property  of  the  definition,  however.  It  is  a  direct  con¬ 
sequence  of  the  far-field  assumption,  whereby  the  effective  coherent 
scattering  area  is  sufficiently  small  vis-a-vis  source  and  receiver  dis¬ 
tances  (RqT,  RqR)  and  dimensions,  cf.  Sec.  5.5  and  Eq.  (5.31)-[4l.  It 
also  depends  on  the  correlation  distances  («-x,Jl  )  of  the  (components  of 
the)  random  wave  surface.  Thus,  is  implicitly  a  function  of  geom¬ 
etry,  where  care  must  be  taken  in  its  use  to  ensure  that  the  conditions 


(A.2-2f) 

(A.2-2g) 

(A.2-2h) 
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governing  the  derivation  of  I ^ ncoh  (=Mx~<X>^^’  (2*7a)-[4],  are  obeyed.* 
A  complete  definition  of  the  incoherent  scattering  cross  section 
includes  the  effects  of  diffraction 


3 


incoh 


E  d(°>  h 

incoh 


+ 


(k) 

incoh  ’ 


(A. 2-3) 


where  the  k  >  1,  are  formally  given  by  (A. 2-1),  with  1^]  . 

i  neon  /  \  f  k )  '  neon 

replaced  by  the  scattering  intensities  I:  ^  cf.  Secs. 


X-<X>' 


7.2C,  8.5  of  [4]. 

Preliminary  estimates  here  of  the  magnitude  of  the  leading  diffrac¬ 
tion  component,  (»  ,  indicate  that  it  is  ordinarily  con¬ 

siderably  smaller  than  the  components  of  the  "classical"  term, 
so  that  in  this  study  we  shall,  at  least  initially,  neglect  the  diffrac¬ 
tion  contributions. 


A. 2-2  Coherent  Scattering  Cross  Sections,  3^ 

The  coherent  scattering  cross  section,  a£j^,  is  formally  the  same 
as  ajpcoh’  (A . 2—1 ) •  except  that  now  the  incident  intensity  (A.2-2b) 
becomes 


rincident(at  the  receiver>coh  *  Ko(0)in/U7T(RoT  +  RoR} 


(A. 2-4) 


which  is  the  "mirror  reflection"  term.  Thus,  I.  .  .  = 

incia-coh 

(1/4)I(ncid-ircoh-  cf-  (A,o2,'2b)-  In  addit1°"'  MncoA  1s  rePlaced  b* 
the  coherent  component  M<^>(0)»  cf.  Sec.  7.2B,  Eq.  (8.7),  [4].  Since  only 

the  "classical"  or  (ka0)  component  of  the  scattered  field, contains  a 


*The  factors  4*  in  (A.2.2a,b),  and  in  (2.9),  (2.10b(  ff . ,  etc., 

arise  because  of  our  definition  of  the  green's  function,  (5.8)-[4],  and 

source  function  (5.3a)-[4]  in  the  equations  of  propagation  (here  a  Helmholtz 

medium).  Thus,  our  acoustic  field  is  =  a/4iT,  where  a  is  derived  from 

a  green's  function  source  of  the  form  -4*6 (R-R' ) 6 ( t- t ' ) .  However, 

-(0) 

because  of  the  particular  form  of  the  definition  of  0  used  here  and 
generally,  the  scaling  of  the  field  is  immaterial,  as  is  required  in 
any  useful  definition. 
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potentially  coherent  contribution  (excluding  any  direct  field  which  may 
be  received  under  certain  mutual  geometries),  the  complete  coherent 
scatter  cross  section  is  now  specifically 


or 


•  /n<  (scattering  at  R)  R2D  ,  , 

Jcoh  5  Iinc1d(at  the  receiver  R)  AR£fr  {path  loss  x  ^9T9R^  } 


(A.2-5a) 


;(o>  .  rqr{4^ror  +  Ror)}2  M<^°) 

COh  (9t9r)2  Aref  KQ(0)in  •  (path  loss) 


(A-2.5bj 


A ,2-3  The  Reference  Area,  AREF 

The  reference  area,  AREF,  appearing  in  the  above  definitions  of  the 
scattering  cross  sections,  (A. 2-1)-,  (A. 2-5),  while  arbitrary,  is  dependent 
on  the  beam  pattern  projections  on  S  .  From  Section  6.5  IV,  (6.56)  of  [4], 
our  choice  of  reference  area  is  specifically 


AREF  =  V2, 


\  =  2*/'AB('S"oT)' 


(A. 2-6) 


where  4^  is  the  projected  area  (on  SQ) ,  here  specifically  of  the  com¬ 
bined  gaussian-omni -directional  beam  patterns  used  in  this  study  and  in 
[4];  A3(0qT)  are  parameters  of  the  projected  beam  pattern  on  SQ,  cf. 
6.3;  also  Secs.  6.4,  6.5  of  [4]  for  details. 


A.2-4  General  Geometry  in  Detail: 

The  Received  Scattered  Wave 

For  the  detailed  development  of  the  analysis  it  is  now  necessary 
to  describe  the  general  geometry  of  the  underice  surface  in  relation  to 
the  source  and  receiver  in  our  current  "one-bounce"  model.  This  geometry 
of  the  general  (far-field)  source-surface-receiver  configuration  is 
shown  in  Fig.  A. 2-2,  in  more  detail  (cf.  Fig.  A. 2-1). 
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P(RT+C ,t) 

I 

'AS  s.underice  surface 


Figure  A. 2-2  Geometry  of  the  transmitter  (@  0T),  underice  surface  (z),  re 
ceiver  (at  0R),  showing  surface  elevation  z  above  ;  =  0:  S  , 
and  the  common  region  of  illumination  (shaded),  on  S  due  0 
to  overlapping  beam  patterns  Ctp  G^;  (here  Vc  =  0).  0 


The  associated  received  scattered  wave  (exclusive  of  diffraction 
terms)  is  (from  (7.7-[4]) 


x(0)'t>  far-field  *  4o  /  S,„(s/2,i )f(0)(5| . . |f0)e“  <£. 


narrow-band 


where  the  Total  Surface  Spreading  Function  (TSSF)  is  now  specifically 


I 


I 


r(0) 


-aw^c^T 
o  o 


RoTRoR  -o 


-( s/c0 ) [  (r+i)*  2Vc0t0]  .  . 

e  dxdy , 


(A.2-7a) 


cf.  (7.7) - [4] ,  subject  to  the  conditions  of  Sec.  5.5,  [4],  summarized 
here  in  A. 2-15  following. 

The  various  elements  of  (A. 2-7)  are  described  below.  These  are: 


A ^  =  effective  "illuminated"  surface  area,  by  joint  beam  pro¬ 
jection  in  S0,  cf.  Sec.  6.3,  Fig.  6.1,  [4]. 

dSQ  =  dr  =  dx'dy'  =  dxdy  on  SQ:  T  =  0;  with  dz  =  dSQ/nz; 

E1  =  x  +  1  ,  where  x'  is  always  on  the  underice  surface  ; 

X  *  i.xx  +  iy;  [ ( x '  ,y‘  )-system  is  simply  a  translation  of  the 
(x,y)-system] 

To  =  (RqT  +  roR)/c0-  time-delay  from  0T  -  0$  -  0R: 

— oT’  *^oR  "-V.oR/ISoT*  ^oR|:  unit  vectors; 

±j’  1r  =  j*y  R/ I Rj  Rl :  unit  vectors,  to  projection  of  scattering 
point  on  Z  upon  SQ; 


(A. 2-8) 

(A.2-8a) 

(A.2-8b) 

(A.2-8c ) 
(A.2-8d ) 
(A.2-8e) 


(A.2-8f) 


RT'  RR’  R0T’"  =  i £rI »  etc.,  distances;  (A.2-8g) 

t '  =  t  -  AtdS  =  doppler-shifted  epoch;  (r^‘  ,t')  are  on  the  under¬ 


ice  surface  z; 


(A.2-8h) 


Atds =  -it  *  i*  (i’t}  =  ij  •  {i +  £(i’t-RR/c0)>  :  rr  4  r0r  -  i0R  •  r;  (A-2_8i } 


L(x,t)  =  (vector)  elevation  of  the  ice-bottom  surface,  vs. 
£  =  0,  the  plane  SQ;  cf.  Figs.  2.2. 


(A.2-8J) 


a  =  t0/2cq  =  absorption  coefficient  (sec  /meter),  cf.  (5.18)-[4]. (A.2-8k) 


Also,  we  have 

CLRj  =  (complex)  beam  pattern  of  R,T;  AQ  =  injected  signal 
ampl itude. 

In  addition,  we  write 


(A. 2-81 ; 
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n2  =  (1  +  ^  cx,  cy  =  §f ,  :  surface  slopes  (A.2-9a) 

_n  =  normal  to  surface  E,  ^Ix^x  +  iy^y  "  j2)nz; 

(inward  normal  =  -_?  )  (A.2-9b) 

v  r>  y.  n  =  iT  of  /c  =  beam-steering  wave  numbers,  cf. 

~oT  oR  T,R  o  o 

(5.13b),  (5.7a),  [4];  (A.2-9c) 

_R_tr  =  Lq  +  lz(hT-hR),  cf.  Fig.  A. 2-2;  vector  distance 

between  Oy  and  0R.  (A.2-9d) 

From  Fig.  A. 2-2,  it  is  readily  shown  that  the  unit  vectors  Jy,  jR 
become  for  these  bistatic  configurations 

jj$  =Jt/|Rt|  =J.x  cos$T  sineT  +  i  sin,j,y  sineT  +  jz  cose  (A.2-10a) 

I  o 

iR|s  =  VI-RrI  =  _{ixRT  C0S*T  Sin6T  +4(RT  Sin<)>T  SineT  "Lo)+  izhR}/RR  * 

(A.2-10b) 

where 

Rr  =  (R2  sin20y  +  L2  -  2RyL0sin<j>y  sin0y  +  h2)\  and  j 

(A.2-10C) 

[F°r  j0j» J0r  we  simply  replace  Ry  by  RQy,  $y  -  iQy.  etc.  in  (A. 2-10).]  We 
also  find  it  convenient  to  write 

'  4<!  +  R(o)T/R{o)R  co“(o)Tsin9(o)T 

*'ly  '0  *  R(o)T/R(o)R>  sl"»(o)TSl"e(o)T  -  VR(o)R! 

+  i2(cosS(o)T  "■  cos0(o)R)-  (A. 2-11) 

cf.  the  exponent  in  (A.2-7a).  [Again,  for  the  reference  vectors  J^y, 

-oR"  ~oT’  etc”  we  set  VR0T’  etc-’  t0  get  2— o  =  ioT  “  -oR * 

as  indicated  in  (A. 2-11).] 
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A  critically  important  simplification  of  our  canonical  result 
(A. 2-7)  results  from  the  ability  to  employ  "narrow  beams"  (.N.B.),  e.g., 
beams  narrow  enough  so  that  the  mutually  "illuminated"  surface  i  (=  S  ) 
(i.e.,  the  shaded  region  in  Fig.  A. 2-2,  for  example)  is  sufficiently  small 


that  the  spatial  geometry  (~Ry,  RR,  etc.)  for  source  and  receiver  change 
little  over  the  "illuminated"  surface  region  (~  S  ).  Thus,  we  can  replace 


— T  by  -oT’  eT  by  0oT 
dependent  quantities 


,  etc.,  and  most  important,  note  that  the  angle- 
in  F^,(A.2-7a)  can  also  be  replaced  by  the  (con¬ 


stant)  reference  quantities  0Qy,  <t>  y,  etc.  This  results  in  a  "factoring" 


of  the  beam  patterns,  in  that  they  now  will  depend  only  on  their  projec¬ 


tions  on  the  surface  about  0^,  and  not  on  the  coordinates  of  the  v 
surface  away  from  0<- .  [For  details,  see  Section  6,  [4].] 


A. 2-5  Parameters  of  Eq.  (3.6):  d:UJ,  R  »1 
- 3 — 1 ‘ 1  nc- — ; 


The  various  elements  of  (3.6)  are  specifically: 


Ro’  Ro  =  mean-square,  mean,  reflection  coefficients  (=1,  for 
water/air  and  water/ice  interfaces); 

!  ,  S  -  mean-square,  mean-,  shadowing  function,  see  Sec.  7.4C- 

|  [4],  and 

2  2  2  2  ^  L 

i  aLx  ’  a  Ly=  <’Lx>’  <^i/> :  (cLx  =  *1T’  etc-):  mean-square  slopes  of  the 

"low-frequency"  ice  surface  component,  cf.  Sec.  7,4£-[4],  and 

(A.2.12J  i  •«$>■■  mean-square  height  (about  <c^>  =  0)  of  the  high- 

frequency  ice  surface  component; 

=  surface  wave  number  intensity  spectrum  of  (cf.  (3.5) 
and  Appendix  (A.l))  with 

55  I  Hf(M0)  cos(k*4r)  » 

G  2ir 

Eq.  (A. 2-6); 

cos  9oT  +  cos  eoR. 
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The  directional  vector  2aQ  is  given  by  (A. 2-11):  for  arbi¬ 

trary  angles  of  illumination  and  observation. 

Important  special  cases  of  (A. 2-11)  are: 

(i) .  Backscattgr  ,C.Rgp:  (1^=  0;  SqR  =  qqT'>  RoR  =  RoT; 

hR  =  hT;  *0R  =  V  +  w/2;  cf.  Fig.  2.2). 

•••  2^  -  2(ix  cos  *qT  sin  9qT  +ly  sin  sin  6qT  +  lz  cos  9qT); 

(A.2-13a) 

(a2  +  ct2  +  a2  )/a  /2  =  2/COS  9  T.  (A.2-13t) 

v  ox  oy  ozy  oz'  °T 

( i i )  .  Bistatic  at  the  Snell  Angle  (R  /  T):  Lq^  0 


^oR^x 

^oT^x 

when  <j>oT  =  n/2  =  UqR  -  n/2);  1 

^  Snel  1 

^>oR^y 

(ioT)y 

when  Lq=  (RqT  +  RqR)  sin  9qT; 

(A.2-14a) 

(-oR}z  = 

'^oT^z 

when  eoR  =  Vt  ] 

|  plane 

2^o  =  : 

’Jz  COS 

eoT’  or  5o  =  iz  cos  9oT’ 

1 

(A.2-14b) 

(aox  +  aoy  +  “ozV^oz72  =  2  cos  90T  *  |  ( A . 2 -  14c ) 


Otherwise,  (A. 2. 11)  is  the  general  relation. 

The  principal  assumptions  and  approximations  pertaining  to  our 
general  high-frequency  results ' (3.6) ,  (3.7),  are: 

Ml).  Far-field  (Fraunhofer)  geometries,  cf.  Sec.  5.5  of  [4],  with  7c=0; 
\(2).  Narrow-band  signals  (so  that  we  may  treat  time  parametrically 
|  in  the  moving  wave  (or  ice)  surface  vis-5-vis  the  acoustic  signal); 
i  cf.  remarks  after  Eq.  (5.21b) -[4]; 

(A.2-15)<^(3) .  Narrow  beams  (cf.  Sec.  6.6-T4];  at  least  one  narrow-beam; 

\(4).  Neglects  diffraction  terms:  k=l:  "Diffuse"  scatter;  k^2: 

!  multiple  scatter  (cf.  Sec.  3.3;  also  Sec.  8.5,  of  [4]); 
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<  (5).  Small  Rayleigh  numbers  for  the  small-scale  surface, 

\  Rh  5  (k0b0aH)2  <<  1; 

\{6).  The  surface  (cH)  is  "small-scale,"  e.g.,  the 

correlation  distance  of  the  low  wave-number  component  is 
,  much  larger  than  that  of  the  high  wave-number  contribution 
(A. 2-15)  <^(7) .  The  small-scale  surface  is  statistically  independent  of 
x  the  c^-surface  component  (vide  remarks  in  Sec.  7.3A,  Eq. 
j  (7.19)  et  seq. ,  [4]); 

(8).  Both  components  of  the  wave  surface  are  essentially  homo- 
j  geneous  and  stationary,  at  least  over  the  "illuminated" 
area  and  for  times  long  compared  to  the  duration  of  the 
incident  signal  (in  the  case  of  platform  motion). 
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Appendix  A. 3:  Second-Order  Probability  Densitigs, 
of  the  Mixed  or  Gauss-Composite  K,S  Process,  Eg.  (2.1) 


Our  task  here  is  to  see  if  it  is  possible  to  generalize  the  mixed  K,S 
process,  described  in  Eq.  (2.1)  by  the  weighted  sum  of  two  first-order 
gaussian  pdf's,  by  a  suitably  weighted  second-order  pdf,  consisting  of  a 
weighted  sum  of  second-order  gaussian  processes. 

We  begin  with  (2.1)  et  seq.,  which  we  rewrite  here  for  the  empirical 
fit  to  the  Lyon  data  for  the  underwater  ice  elevation  c,  as  the  first- 
order  pdf 


Wj(5)  - 


bj  -  (z;-ms)2/2o| 
^2t to  2 


b2  -U-mK)2/2a£ 
/2itok 


(A. 3-1) 


for  which  specifically 

b1  =  0.930;  b2  =  l-bx  =  0.070;  (b^  =  1)  \ 

m$  =  4. 50(m) ;  m^  =  I2.5(m)  l  (A. 3-la) 

=  l.OO(rn);  a ^  =  3.00(m) 

The  associated  characteristic  function  (c.f.)  (2.4),  or  (A. 3-1)  is  therefore 


F1(i^)  =  bje 


imsC-a^2/2  im^-aj^2^ 
+  b2e 


=  b  p(  1 )  +  b  F^ 

r  IS  VlK 


(A. 3-2) 
(A. 3-2a ) 


Accordingly,  we  obtain 
:( 1 ) 


dF 


C  =  -i 


1 


d£ 


5=0 


=  b1m$ 


T.  d‘Fl 

5  ‘ 


=  b1(os 


5=0 


+  b2mK  =5.1  (m); 

+  m2)  +  b2(o2  +  m2)  =  31.3  (m2); 


(A.3-4a) 


^  *  .  *  .  *  .  ,  **  a 


(A.3-4b) 
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•*.cr2  =  £2  -  =  var  c  =  b^2  +  b,^2  +  bjb2(ms  -  m^)2  =  (2.38)2  (m2), 

(A.3-4c) 

since  b^  2  =  1  "  b2  1*  cf*  (A-3-la)- 

In  our  effort  to  construct  a  second-order  pdf,  w2(c  1  »tj  ’ 

which  is  a  valid  generalization  of  the  first-order  pdf  (A. 3-1),  we  must 
require  that  the  following  conditions  be  obeyed: 

I.  The  first  and  second  moments  derived  from  w 2  (or  the  associated  c.f., 
p£_^)  are  the  same  as  those  obtained  from  the  first-order  w^, 

(A. 3-1),  viz.,  (A.3-4a,b),  with  the  same  variance  (A.3-4c),  e.g., 
var  =  var  C2  =  var  5 ,  for  C  now  postulated  to  be  homogeneous. 

[if  this  latter  condition  is  removed,  then  the  form  of  var  2  is 
to  be  the  same  as  that  of  var  5.] 


II.  When  Ar  =  ^2=-Ki  (1n  the  (x»y)~P^ane  of  t(r))  becomes  infinite,  i.e., 
A r  -*■  +<*>,  then 


lim  F  =  f(1)f(2) 
Ar-*±°°  2-C  1-C1-C 


Fi(iWi(u' 


(A. 3-5) 


c2  are  (second-order)  independent  asAr  +  +<»,  where  F^jr),  is 
given  by  (A. 3-2)1 


Condition  II  then  requires  that  the  second-order  c.f. 


form 


F,  have  the 


(A. 3-6) 


where  F.<-  K  has  the  form  (A.3-2a),  for  inasmuch 
f1s!ikf1s!ik'  "e  see  that 


as 


lim  F 

Ar->±®  2S,2K 


lim  F  =  (b  F(1) 
Ar-*±°°  r2-c  1  r  IS 


+  b  r(l)wb  r(2) 
b2FlK  }(blFlS 


+  b  p^2^ )  =  F^M2^ 

V2K  1  U-C  1-t 


(A. 3-7) 


with  f|^,  etc.,  given  formally  by  (A.3-2a).  Note  that  no  specific  stat¬ 
istical  structure  for  Us,cK,0  has  yet  been  imposed.  [Also  note  that 

^2-^1^  =£2=0  =  ^bl  +  2blb2  +  b2^  =  ^bl  +  b2^  =  lr*  as  re9uired-] 
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Now,  let  us  require  that  be  given  specifically  by  the  gaussian 
first-order  composite  form  (A. 3-2)  and  that  c 5 .  and-*-£,  are  homogeneous 

processes.  Then  the  (first-order)  first-  and  second  moments  are  given  by 
(A. 3-4),  as  required.  However,  according  to  (A. 3-7)  this  has  been  es¬ 
tablished  only  for  in  the  limit  Ar We  need  now  to  determine 
whether  I  (and  II)  are  satisfied  in  the  general  case  0  <  |at|  <  ±®. 
Accordingly,  for  ^2-^  to  reduce  to  the  gaussian  form  (A. 3-2)  we  require 
that  be  second-order  gaussian  processes,  so  that  (A. 3-6)  becomes 

explicitly  (for  these  homogeneous  cases) 


2  imsf’l+imS^2"°S^l+^2+2pS^1^2^ 

F2-c  =  ble 

+  b  b  e1mS^l+imK52"^aS?l+aK52+0SaKpSK51^2+aSaKpKS^1^2^ 


+  blb2e 


1mK^1+imS^2"is[aK^l+aS^2+0S0KPKS^1^2+0SaKpSK^1^2^ 


2  imK5l+lmKC2"J50K(ClH2+2pK^1^2) 

’’’  02=  • 


(A. 3-8) 


Tim 

Clearly,  since  m  Ps’<3ks,pSK~>’®’  (^.3-8)  reduces  directly  to  the  form 
(A. 3-7),  as  required  by  II,  (A. 3-5).  Moreover,  -1 «£ Ps»pks,pSK^  these 
normalized  auto-  and  cross-covariances  are,  by  definition,  bounded  as  in¬ 
dicated.  In  addition,  the  quadratic  forms  in  the  various  exponents  of 
(A. 3-8)  must  be  at  least  positive  semi-definite.  This  means  that  no 
adjustable  scale  factors  can  be  introduced  as  coefficients  of  p^p^.  etc., 
or  of  a5>o2,  etc. 

Next,  let  us  obtain  the  various  second-order  moments  associated  with 

(A. 3-8).  Inasmuch  as  F,  here  reduces  to  (A. 3-7),  with  f|^,  etc.,  given 

7  2  ^  *  £ 

by  (A. 3-2),  s  ^2 =  5  ’  (A.3-4b),  in  these  homogeneous  cases.  For  the 

covariance  K^.(Ar)  we  get  directly  on  differentiation  of  (A. 3-8), 


Kc(£n)  =  ?x?2  '  W  = 


34 


2-; 


Ci=C2=0 


-  (b1ms+b2mK)2  (A.3-9a) 


=  b2a^ps  +  b2a2P)<  +  b^o^p^+p^) , 


(A.3-9b) 
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r- 


which  vanishes,  as  required,  when  Ar-*-±°°,  since  ,  Pc,  (etc.)-*-0. 

Ar  -♦'Zoo  o 

We  can  rewrite  (A.3-9b),  using  =  l-b2>  etc.,  as 

Kr-(  Ar )  =  b ^ ^  b2°KpK  +  b  ^  b2 ^PS ^  PSK^ PKS ^  ^  *  (A. 3-10) 

The  behavior  as  Ar-*0  now  needs  to  be  examined.  We  have  for  (A. 3-10) 


iy°)  ■  b/s  *  b2°K  +  blV2°sV,>  -  4  ■  4)  -  c2  -  c2; 


P  =  „(0)  ,  JO) 
Po  "  PSK  PKS  • 


(A. 3-11) 


This  is  to  be  compared  with  C  calculated  directly  from  the  first-order 
c.f.  (A. o-2),  viz.  (A.3-4c) : 


(var  c  =  42  -  C2) 


lst-order 


bl°S  +  b2°K  +  b1b2(ms_mK)2-  (A. 3-12) 


We  see  at  once  that  it  is  quite  impossible  for  (A. 3-11)  and  (A. 3-12)  to  be 
equal  (bj  f  b2  f  1.0)  unless  m$  =  and  p^  =  pj^  =  +1,  with  aK  =  o^: 
then  the  two  (K,S)  processes  are  second-moment  identical,  and  correlated. 

This  is  not  at  all  physically  likely,  so  we  must  discard  it  as  an  accept¬ 
able  model . 

We  note  that  changing  the  weighting,  b^b2  (in  such  a  way  as  to  obey 
II,  cf.  (A.3-7)),  does  not  help  us  achieve  agreement  between  var  ^st 
(A. 3-12)  and  K^(0),  (A. 3-11).  For  example,  let  us  replace  b^b2  by 
b1b2(l-pSK),  b1b2 ( 1- )  respectively  in  (A. 3-8).  This  gives  the  correct 
form  for  Ar-*-±°°,  (A.3-7),  but  merely  inserts  a  factor  l-P^.  1-PK$  with 
PSK’  PKS  in  (A’3-10),  (A. 3-11),  which  again  does  not  allow  an  equality 
between  (A. 3-11)  and  (A. 3-12),  except  in  the  trivial  case  noted  above. 

However,  if  we  set  b^=l  (-'-b^O),  or  vice  versa,  we  do  get  agreement, 
as  we  would  expect,  since  only  one  gauss  component  is  present.  Of  course, 
the  choice  b^=l  (b2=0)  gives  a  good  fit  to  the  first-order  data  except 
for  the  important  tails  (represented  by  the  K-gauss  process  here).  In 
retaining  this  model  (Model  I),  we  make  the  comparison  indicated  in  Sec. 

2.1,  Eq.  2.7:  we  choose  here  a  single  gauss  process,  4,  whose  mean  and 
variance  are  those  of  the  gauss  composite  process  obeying  (2.1),  or  (A. 3-1)- 
(A.3-4)  above.  Thus  our  Models  Ia,b  are  taken  to  be  a  single  process,  with 
first-order,  first-  and  second-moments  determined  from  the  Lyon  data  fit  (2.1) 
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Before  we  note  the  implication  of  the  above,  we  observe  that  using  a 
composite  gauss  process 

pcomp  =  Vs  +  a2;K’  (0  <  alfa2  <  1),  (A. 3-13) 

rather  than  the  gauss  composite  process  (described  by  (A. 3-1),  (A. 3-2)), 
also  does  not  permit  an  equivalence  between  the  two  (unless  a2  =  0,  for 
example).  This  is  easily  seen,  for 


"comp  =  V’S  *  V*’  "  V  ;S  =  mS;  d2  =  b2;  ;K  =  mK 

(0  d  brb2  <  1)  (A. 3-14) 

Similarly,  we  get  dmect!y 

var  ’comp  =  Vs  *  Vr.  *  blb2  ,  2KSK(0)  ‘  JS  ‘  3K};  (A. 3-15) 

ksk(0)  »  msk(0)  -  m^m^.  MSj<(0)  1  ;$4K-  (A . 3- 15a ) 

From  (A.3-4c)  and  (A. 3-15),  we  can  write 

var  ’comp  =  °l  *  blb2  {oK  +  °S  '  2oS°KpKS(0)  +  (nVmK)2}  >  °’  (A'3'16) 

with  pks(0)  =  PSK(0)  =  ksk(°)/0sjk'  Consequently, 

var  ^comp  *  °2  >  blb2  lo<  *  °s  -  2oS°K0KS(0>  *  ("s-”k>2)  ’  °-  (A-3'17> 

2 

so  that  var  Ccomp > (>  0)  always:  these  two  first-order  processes  (com¬ 
posite  gauss,  (A. 3-13),  and  gauss  composite  (2.1))  are  not  the  same.  We 
need  proceed  no  further  in  the  second-order  cases:  the  second-order  forms 
cannot  reduce  to  a  common  first-order,  (2.1). 

In  summary,  our  principal  results  here  and  their  implications  are 
concisely: 

I.  It  is  not  possible  to  construct  a  second-order  gauss  composite  process 
which  reduces  to  the  prescribed  first-order  process  [(2.1)],  having 
the  same  second-order  (and  therefore,  higher-order)  moments. 
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II.  It  is  possible  to  construct  a  single-component  process  with  the  same 
first-  and  second-  (first-order)  moments.  This  is  the  procedure 
followed  in  the  substance  of  this  Report,  for  which  Model  I  is  thus 
essentially  treated  as  a  single-component  process,  cf.  (2.16)  ff . , 
and  (2.7)  et  seq.  The  moments  used  explicitly  are  obtained  by  a 
combination  of  fitting  to  the  data  and  physical -model  building. 
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